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I  INTRODUCTION 


1.1.  Importance  and  Problem  Definition 

Precision  frequency  sources  such  as  quartz-crystal  oscillators  and 
atomic  frequency  standards  and  clocks  find  a  wide  variety  of 
applications  in  many  engineering  fields  [1.1]  ,[1.2].  Some  important 
application  areas  are  doppler  radar  systems;  satellite  and  spacecraft 
for  guidance,  tracking  and  communications;  range  measurements  and 
digital  communication  systems.  In  these  and  other  coherent  information 
processing  systems,  stability  of  the  precision  frequency  sources  is  of 
prime  importance  in  maintaining  synchronization.  Most  oractical 
sources,  however,  exhibit  some  degree  of  frequency  instat  ty.  This 
may  be  due  to  internal  or  external  environmental  conditions  which  the 

frequency  source  is  exposed  such  as  additive  and  multiplic  .a  noise, 
temperature  fluctuations,  supply  voltage  variations,  magnetic  field, 
atmospheric  pressure,  humidity  and  vibration.  These  and  other  time 
dependent  physical  causes  generate  random  time-varying  fluctuations  and 
instabilities  in  the  parameters  which  characterize  the  oscillator,  viz., 
frequency  and  phase,  which  result  in  a  degradation  in  the  performance  of 
the  coherent  information  processing  systems.  In  order  to  evaluate  the 
actual  system  performance,  an  accurate  characterization  of  the 
instability  and  its  description  in  terms  of  some  quantitative  measures 
is  essential  [1.3],  [1.4].  By  phase /frequency  instability,  we  mean  the 
degree  to  which  an  oscillator  produces  random  unwanted  and  time 
dependent  departures  from  a  pure  sine  wave.  In  the  spectrum 

representation,  it  is  the  degree  to  which  these  departures  cause  the 
oscillator  specular  component  to  disperse.  The  consequences  of 
instability  are  losses  in  temporal  synchronization  and  spectral 
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resolution  in  the  associated  system.  In  this  dissertation,  we  generally 
model  the  phase  error  process  as  power  law  spectral  density  models 

[1.3] ,  i.e.,  power  spectral  density  is  modeled  as  a  polynomial  in  the 
frequency  domain.  These  models  tend  to  have  certain  physical  origins, 

[ 1.4]  ,  as  follows: 

a)  White  phase  noise,  y(u),  is  usually  due  to  white  noise  sources 
external  to  the  oscillator  loop. 

b)  Random  walk  phase  (white  frequency)  noise,  x(u),  usually 
arises  from  additive  white  noise  internal  to  the  oscillator 
loop . 

c)  Random  walk  frequency  noise,  v(u),  is  usually  related  to  the 
oscillator  environment  such  as  temperature,  vibration  and 
shocks . 

The  objective  here  is  to  characterize  these  undesired  time-dependent 
frequency  departures  and  then  utilize  it  for  the  evaluation  of  their 
effects  on  the  performance  of  coherent  information  processing  systems. 

Since  any  frequency  variation  implies  an  associated  phase 
variation,  we  shall  concentrate  on  the  phase  random  process  as  the 
source  of  frequency  instability  and  develop  stochastic  models  for  it. 
Although  the  model  and  the  related  discussion  is  valid  for  numerous 
applications ,  emphasis  will  be  placed  on  radar /communication 
applications.  In  these  applications,  frequency  instability  is  just  one 
of  the  many  sources  of  system  performance  degradation.  Frequency 
instability  induces  a  spectral  spread  which,  in  turn,  affects  the 
probability  of  detection  of  the  radar /communication  system. 


This 


dissertation  treats  only  the  effects  of  phase /frequency  instability  on 
the  system  performance.  Other  sources  of  system  performance  degradation 
such  as  additive  and  multiplicative  noise;  doppler  spread  of  desired  and 
undesired  scatterers;  atmospheric  chatacteristics;  spurious  modulations 
and  nonlinearities  from  other  components;  channel  mismatch  and 
fluctuations  etc.  are  not  discussed  here.  They  have  been  discussed  in 
detail  elsewhere,  e.g.  [1.5]— [1.8] . 

1.2 .  Previous  Work 

As  Indicated  above,  characterization  of  frequency  instability  is  an 
important  problem  and  has  received  considerable  attention.  This  problem 
was  addressed  extensively  in  a  special  issue  of  the  Proceedings  of  the 
IEEE  [1.2].  Throughout  that  issue,  frequency  instability 
characterization  in  both  frequency  and  time  domains  and  translations 
between  them  were  considered.  Topics  in  frequency  control,  frequency 
stability  requirements  in  various  fields  of  applications,  and  progress 
in  atomic  and  quartz  frequency  control  were  discussed.  That  issue  was 
an  outgrowth  of  the  tremendous  effort  put  forth  in  the  United  States 
space  program  to  satisfy  the  urgent  need  to  standardize  terminology.  It 
was  a  first  step  toward  bringing  together,  in  a  coordinated  way,  the 
latest  information  related  to  this  field  including  a  number  of  papers 
dealing  with  the  fundamentals.  There  still  remained  individual 
preferences  for  terminology  and  notation.  Three  classes  of  papers 
constituted  the  issue:  papers  dealing  with  theroetical  issues,  papers 
providing  the  details  of  oscillators,  and  papers  describing  the 
requirements  and  specifications  from  the  user's  point  of  view. 


Cutler  and  Searle  [1*9]  discussed  sources  of  quartz  oscillator 


noise  and  associated  frequency  fluctuation  with  respect  to  averaging 
time.  Atomic  standards,  both  passive  and  active,  using  servo-controlled 
oscillation,  were  discussed  with  respect  to  servo  time  constants. 
Conventional  methods  for  handling  random  signals  were  applied  to  phase 
and  frequency  fluctuations  for  obtaining  meaningful  criteria  for 
specifying  oscillator  stability. 

Vessot  et  al  [1.10]  described  an  auto-correlation  technique  for 
measuring  the  short-term  properties  of  stable  oscillators.  The  phase 
and  frequency  were  converted  to  amplitude  fluctuations  via  a  function 
multiplier,  a  second  oscillator,  and  a  time  averager.  Mean  square 
averages  of  frequency  phase,  and  amplitude  were  obtained  from 
measuremants  on  hydrogen  masers. 

Allan  [1.11]  developed  the  theoretical  relationship  for 
determination  of  frequency  standard  deviation  for  any  finite  number  of 
data  samples  with  respect  to  the  infinite  time  average.  He  also 
provided  a  method  for  determining  power  spectral  density  which  is  system 
dependent,  i.e.  it  depends  on  bandwidth  and  sampling  time.  Such 
diverse  topics  as  flicker  noise,  masers,  the  cesium  beam  standard,  and 
rubidium  gas  cells  were  also  discussed.  These  topics  were  further 
addressed  in  [1.3].  A  measure  of  frequency  stability  was  defined  with 
respect  to  the  spectral  density  of  fractional  frequency  fluctuations. 
An  alternate  measure  was  also  given;  the  expected  variance  for  the  N 
sample  averages  of  the  fluctuations  taken  over  some  time  duration,  t, 
with  the  begining  of  successive  periods  spaced  every  T  units  of  time. 
The  variance  was  denoted  by 

2 

<r  (N,T,t) 

y 
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A  preferred  measure  of  frequency  stability  in  the  time  domain  was 
achieved  by  choosing  N  ■  2  and  T  ■  t  to  obtain  the  "Allan  variance", 
which  found  widespread  use  in  situations  where  system  noise  bandwidth 
enters  into  the  measurements.  Systematic  variations,  frequency  and  time 
domain  relationships,  applications,  and  measurement  techniques  were  also 
treated. 

An  excellent  survey  of  the  state-of-the-art  in  this  area  was 
presented  in  [1.4].  Many  parameters  proposed  for  phase  and  frequency 
instability  characterization  were  reviewed  including  both  the  widely 
used  classical  concepts  and  the  more  recent  less  familiar  approaches. 
They  included  those  recommended  by  the  IEEE  Subcommittee  on  Frequency 
Stability  and  by  Study  Group  7  on  "Standard  Frequencies  and  Time 
Signals"  of  the  International  Radio  Consultative  Committee  (CCIR).  The 
paper  treated  transformations  of  the  error  phase  and  frequency  functions 
between  time  and  spectral  domains  and  the  advantages  of  new  approaches 
that  provide  improved  understanding.  Here,  the  structure  functions  of 
phase  and  its  relations  with  the  sample  variance  and  the  Hadamard 
variance  were  discussed.  Polynomial  phase  and  frequency  drifts  along 
with  random  processes  were  modeled  by  power-law  spectral  densities. 

Specific  frequency  stability  issues  related  to  radar  have  been 
presented  by  Raven  [1.12]  who  derived  oscillator  requirements  for 
coherent  radars.  The  modulation  effects  imposed  by  phase  instabilities 
within  the  oscillator  were  examined  for  their  effects  on  the  radar 
waveform  and  system  performance.  Leeson  and  Johnson  [1.13]  treated 
short-term  frequency  stability  along  with  system  and  circuit 
requirements  for  operation  of  Doppler  radars  in  severe  vibration  and 
acoustic  environments.  Linewidth  and  spectrum  are  used  to  define  the 
short-term  stability  of  a  Doppler  radar. 


Lindsey  and  Meyr  [1.14]  determined  the  time  dependent  probability 
density  function  of  the  phase  error  as  well  as  the  distribution  of  cycle 
clips  using  a  Wiener  (Brownian  Motion)  process  embedded  in  a  renewal 
process  of  a  first  order  correlative  tracking  system  with  periodic 
nonlinearity.  The  relation  between  phase  and  instantaneous  frequency 
was  studied  by  Boileau  and  Picinbono  [1.15].  In  particular,  conditions 
under  which  the  phase  cannot  be  stationary  or  the  instantaneous 
frequency  does  not  exist  were  discussed.  A  new  method  of  analysis,  the 
method  of  finite-time  frequency  control  to  describe  the  frequency 
stability  was  compared  with  the  Allan-Variance  procedure. 

The  Allan  variance  was  used  by  Yoshimura  [1.16]  to  estimate  the 
variance  of  several  simulated  phase  random  processes  containing 
frequency  noise.  Distributions  of  the  estimates  were  compared  to  the 
chi-square  distribution.  It  was  found  that  the  auto-correlation 
function  of  the  process  affects  the  outcome  for  white,  flicker,  and 
random  walk  phase  noise.  For  flicker  and  random  walk  frequency,  the 
process  can  be  regarded  as  nearly  independent  for  a  reasonable  estimate 
of  the  variance. 

In  several  papers,  Lindsey  and  Lewis  [1.17]  and  Lindsey  and  Chie 
[1.18]  used  structure  functions  in  characterizing  and  relating  phase 
Instability  with  time  and  frequency  domain  measures.  They  demonstrated 
that  the  fractional  frequency  deviation  Introduced  by  Cutler  and  Searle 
is  related  to  the  first  order  phase  structure  function.  Higher  order 
structure  functions  in  phase  and  frequency  were  also  treated  for  their 
usefulness  in  overcomming  problems  associated  with  frequency  drift  and 
flicker  type  noise  convergence. 


In  this  dissertation,  a  covariance  matrix  approach  for  the 
characterization  of  phase  instability  is  developed.  This  allows  a 
uniform  and  systematic  treatment  of  all  the  above  models  and  their 
properties  which  include  those  of  Lindsey's  and  Chie's  structure 
functions.  The  approach  provides  a  geometric  mathematical  framework 
which  allows  for  the  concise  establishment  and  organization  of  the 
stochastic  operations  required  for  the  treatment  of  auto-correlation 
functions,  characteristic  functions,  and  power  spectral  densities  of 
oscillator  rf  signals  having  phase  instabilities.  We  identify  a  class 
of  non-stationary  phase  random  processes  which  allow  determination  of 
the  rf  signal  power  spectral  density.  We  also  show  that  under  certain 
conditions  on  these  processes  the  mean,  auto-correlation  function,  and 
power  spectral  density  of  the  rf  signal  can  be  determined.  Stationarity 
and  ergodicity  conditions  for  the  rf  signal  are  given  through  a  number 
of  theorems.  Also,  since  we  are  interested  in  communications  and  radar 
applications,  we  consider  the  effects  which  the  oscillator  phase 
instabilities  have  on  modulation  with  envelope  waveforms.  The  total 
effects  due  to  phase  instability  and  envelope  modulation  are  illustrated 
through  numerous  examples. 

1.3.  Analytical  Background  and  Approach 

To  determine  the  spectral  spreading  of  a  radar  or  communications 
waveform  as  affected  by  oscillator  phase  instabilities  we  call  to  mind  a 
number  of  analytical  relationships.  We  make  use  of  the  fact  that  the 
spectral  characterization  involves  the  Fourier  transform  of  the 
auto-correlation  function  or  average  auto-correlation  function,  Rf(u), 
of  the  waveform,  f(u).  Such  a  waveform  is  the  product  of  the  system's 


deterministic  modulating  envelope  waveform,  fo(u),  and  the  oscillator  rf 

signal,  b(u),  which  in  general  contains  the  stationary  and 

non-stationary  phase  instabilities,  z(u).  It  will  be  shown  later  that 

b(u)  is  stationary  in  auto-correlation  function  for  a  class  of 

non-stationary  random  processes,  z(u).  By  the  oscillator  rf  signal, 

_1 

b(u),  we  mean  [1.17]  the  sinusoidal  function  of  the  phase  random 
process,  say  z(u),  which  is  a  function  of  time,  u.  Hence 

b(  u  )  -  A  exp  {  j  [0  +  z(  u  )]  }  (1*1) 

where  A  and  are  arbitrary  constants.  Modeling  of  b(u)  is  considered 
in  detail  in  chapter  two.  The  random  process,  b(u),  may  contain  the 
effects  of  short,  intermediate  and  long  term  oscillator  phase 
instabilities  as  well  as  other  non-deterministic  modulations  induced  by 
propagation  effects,  thermal  drift,  and  doppler  spreading  of  scatterers. 
In  this  dissertation,  we  will  concentrate  only  on  the  effects  of 
oscillator  phase  instabilities.  System  average  power/energy  spectrum 
will  be  determined  for  a  variety  of  waveforms  affected  by  a  set  of  phase 
instability  parameters  using  the  procedures  as  described  below. 

Let  us  denote  by  f  (u),  a  sample  function  of  the  modulated  signal 
k  •  2 

obtained  when  we  multiply  the  modulating  envelope  waveform  ,  f  (u), 

o 


The  function  indicated  in  equation  (1.1)  is  also  known  as  the 
cisoidal  function.  In  this  dissertation,  the  term,  "sinusoidal",  is 
used  to  describe  this  function. 

2. 

The  tilda  (“)  over  the  function  denotes  complex  low  pass 
equivalent.  For  notational  simplicity,  this  symbol  will  be  omitted  with 
the  assumption  that  the  real  nature  of  the  phase  random  process,  z(u), 
and  the  complex  nature  of  the  rf  waveforms,  f(u)  and  b(u),  is  implied 
throughout  this  dissertation. 
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(which  in  general  is  non-stationary )  by  a  sample  function,  b  ( u ),  (which 

k 

is  a  stationary  random  process,  b(u))  i.e., 


f  (u)  =  f  (u)  b  (u) 
k  ok 


(1.2) 


Figure  1.1  shows  the  equivalent  block  diagram  of  such  a  modulation 

operation.  It  has  been  shown  [1.19],  [1.20]  that  the  time 

auto-correlation  function  of  this  product  may  be  determined  as  follows. 

First  we  write  the  time  average  auto-correlation  function,  R  (  t ),  of 

fk 

the  modulated  signal.  For  a  power  limited  waveform,  we  write 


R  (t) 

fk 


lira 

T — »ap 


-  (  £<u 

2T  J  k 


)  f  (  U-t  )  du 
k 


lim 
T — >CB3 


i 

1  /  *  * 

\  b  (u)  b  (u-t)  f  (u)  f  (u-t) 
2T  /  k  k  o  o 


du  (1.3a). 


and  for  an  energy  limited  waveform  we  write  , 


R  (t) 

fk 


oo 

l  f  ( u  )  f  ( u-t )  du 
)  k  k 


b  (u)  b  (u-t)  f  (u)  f  (u-t)  du 
k  k  oo 


(  1.3b). 


We  will  encounter  both  mean  square  power  limited  and  mean  square 
energy  limited  modulated  waveforms.  Without  loss  in  generality  we  will, 
in  this  chapter,  adhere  to  the  notation  for  the  mean  square  energy  limited 
case  thereby  omitting  division  by  2T  and  the  associated  limiting  operation. 
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FIGURE  1.1.  MODULATOR  BLOCK  DIAGRAM 
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The  Foufier  transform  of  R  (  t  )  in  equation  (1.3a)  and  equation  (1.3b) 

fk 

is  the  power  spectrum  and  energy  spectrum  of  the  waveform  respectively. 

These  functions  are  complex  lu*  pass  equivalents  and  *  denotes  complex 

conjugate.  Because  the  time  domain  process,  f  (u),  and  likewise,  f  (u) 

o  k 

defined  in  equation  (1.2),  is  generally  not  stationary,  R  (t)  defined 

fk 

in  equation  (1.3)  is  not,  in  general,  equal  to  the  ensemble  average 
auto-correlation  function,  i.e. 


R  (t) V-  E  {  f  (u)  f  (u-t)  } 
fk  k  k 


Therefore,  R  (  t )  is  not  the  inverse  transform  of  a  power  spectral 
fk 

density.  In  fact,  R  (t)  obtained  through  equation  (1.3)  will  be 

fk 

different  for  different  sample  functions,  b  (u),  and  thus  it  is  a  sample 

k 

function  of  a  random  process  itself.  Since  R  ( t )  is  a  sample  function 

fk 

from  a  random  process,  we  can  determine  the  expected  value  of  R  ( t  ), 

fk 

(likewise  the  Fourier  transform,  S  (w))  by  taking  the  ensemble  average. 

fk 


R  (  t  )  -  E  {  R  (  t  )} 
f  k  fk 


*  it 

b  (u)  b  (u-t)  f  (u)  f  (u-t)  du} 
k  k  o  o 


-a© 


Interchanging  the  order  of  averaging  and  integration  we  obtain 


OP 


*  * 

R(t)  =  \  E{b  (u)  b  (u-t)}  f  (u)  f  (u-t)  du 


k  k  k 


o  o 


-o® 


c© 


S' 


(t)  f  (u)  f  ( u-t  )  du 
boo 


-ao 


The  result  is  the  average  auto-correlation  function  of  f(u)  [1.20], 
section  12.3. 

(1.4) 


R  (t)  -  R  (t)  R  (t) 
f  bo 


QtD 


where 


( 

f(u)f  (u-t)  du  is  the  time  average 
o  )  o  o 


-a© 


auto-correlation  function  of  the  modulating  envelope  waveform,  f  (u), 

o 

and  R  (t)  is  the  ensemble  average  auto-correlation  function  of  the 
b 

process,  b(u).  Therefore,  the  expected  value  of  the  time  average  sample 


auto-correlation  function,  R  ( t ),  i.e.  the  average  auto-correlation 

fk 

function  is  the  product  of  the  stationary  process  auto-correlation 


function  and  the  time  average  envelope  auto-correlation  function.  Upon 


modulation,  the  carrier  shifted  spectrum  of  f  (u)  is  contaminated  by 

o 

undesired  multiplicative  factors  which  are  contained  in  the  oscillator 


signal,  b(u).  The  average  power /energy  spectrum  of  the  product  of  the 


envelope  waveform,  f  (u),  and  the  random  process,  b(u),  is  the  Fourier 

o 

transform  of  the  average  auto-correlation  function,  R  ( t ),  given  in 

f 

eq*.._tion  (  1.4  ). 


The  Fourier  transform  of  the  average  auto-correlation  function  is 
denoted  as  the  average  power  spectrum  o if  the  average  energy  spectrum 
depending  on  whether  the  waveform  is  power  limited  of  energy  limited, 
respectively.  Further  discussion  can  be  found  in  [1.20]. 
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1 


< 


1 
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The  auto-correlation  function,  R  ( t  ),  which  we  have  described  in 

f 

equation  (1.4)  can  be  represented  as  a  product  of  auto-correlation 

functions  corresponding  to  other  independent  random  processes  [1.21], 

b(u),  b(u), .  b(u);  each  representing  independent  phase 

12  n 

instability  processes.  The  average  auto-correlation  function  will 
then  be  [1.20] ,  [1.21] ,  [1.22] 

n 

R  (t)  -  R  (t)Tf  R  (t)  (1.5) 

f  o  lib 

n  i-1  i 


1.4.  Organization  of  the  Dissertation 

In  chapter  two,  we  consider  modeling  and  characterization  of  the 
phase  error  process  and  its  relationship  to  the  frequency  stability 
measures.  Models  developed  include: 

i)  a  superposition  of  a  stationary  white  and  a  non-stationary 
random  walk  phase  process  and; 

ii )  a  random  walk  frequency  process  along  with  a  random  drift 
phase . 


The  white  phase  and  random  walk  phase  models  will  be  treated  together  in 
one  model  while  the  random  walk  frequency  will  be  treated  separately  in 
another  model.  Thus,  the  first  model  treats  the  effects  of  white  noise 
on  the  oscillator  while  the  latter  characterizes  the  oscillator 
instability  caused  by  the  oscillator  environment  [1.4].  It  is  shown  how 
moments  of  the  rf  signal  random  process  including  the  auto-correlation 
function  can  be  evaluated  through  the  use  of  the  phase  covariance  matrix 
and  the  characteristic  function. 
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In  chapter  three,  conditions  are  given  for  stationarity  and 
ergodicity  of  the  mean,  auto-correlation  function,  and  power  spectral 
density  of  the  rf  signal.  Each  of  the  phase  instability  models 
mentioned  earlier  are  treated  here.  Wide  sense  stationarity  is 
established  by  providing  conditions  under  which  the  mean  and 
auto-correlation  function  of  the  rf  signal  are  invariant  over  the  time 
variable,  u,  for  those  phase  instability  models  considered.  A  number  of 
theorems  from  Papoulis  [1.20]  pp  328,  330,  343-344,  are  utilized  to 
determine  the  ergodicity  of  our  models.  It  is  found  that  the  statistics 
of  the  phase  appearing  the  instant  at  which  the  observation  commences 
plays  an  important  role  in  determining  the  stationarity  and  ergodicity 
of  the  rf  signal.  For  Gaussian  distributed  instabilities,  i.e.  white 
phase,  random  walk  phase,  and  random  walk  frequency  processes, 
conditions  at  the  time  the  rf  signal  is  first  observed  are  determined  in 
establishing  ergodicity. 

Spectral  spreading  due  to  the  oscillator  instabilities  in  an 
envelope  modulated  signal  is  discussed  in  chapter  four.  Three  important 
modulating  waveforms  are  considered:  CW,  Infinite  Pulse  Train,  and 
Finite  Pulse  Train.  Spectral  properties  are  determined  by  taking  the 
Fourier  transform  of  products  of  the  auto-correlation  functions 
developed  in  chapter  two  and  time  average  auto-correlation  functions  of 
the  particular  modulating  waveform.  In  the  case  of  random  walk 
frequency  instability,  asymptotic  expansions  are  developed  to  evaluate 
the  related  Airy  integral  which  arises  in  the  determination  of  the 
spectral  spreading. 

In  chapter  five,  numerical  results  and  some  approximations  are 
presented.  An  example  for  each  of  the  major  phase  models  and  modulating 


envelopes  is  worked  out  and  the  results  are  plotted.  Finally,  our 
results  are  summarized  along  with  the  significant  conclusions  in  chapter 
six. 

The  details  of  many  of  the  derivations  are  included  in  four 
appendices.  In  appendix  A,  it  is  shown  how  to  determine  the  n  th  order 
moments  for  the  rf  signal,  b(u).  The  procedure  is  based  on  finding  the 
characteristic  function  of  the  vector  phase  random  variable,  z,  and 
setting  its  argument  to  the  appropriate  value. 

In  Appendix  B,  we  determine  the  covariance  matrix  of  the  phase 
random  process  consisting  of 

1)  the  white  phase  random  process  and 

2  )  the  random  walk  phase  random  process 

together  as  a  superposition.  The  results  of  Appendix  B  are  used  to 
determine  oscillator  auto-correlation  functions  in  chapter  two. 

Appendix  C  is  devoted  entirely  to  finding  the  covariance  matrix  of 
the  phase  random  process  resulting  from  a  random  walk,  i.e.  Wiener, 
frequency  process.  Due  to  the  mathematical  difficulties  encountered  in 
undertaking  such  a  task,  three  approaches  were  used;  the  approximating 
function,  the  integration  formula,  and  a  state  space  approach.  The 
results  of  all  three  approaches  were  in  agreement.  Again,  these  are 
used  in  chapter  two  to  determine  auto-correlation  functions. 

In  order  to  prove  the  theorems  stated  in  chapter  three,  higher 
order  moments  of  the  rf  signal  process  had  to  be  determined  for  the 
phase  instability  models  considered  here.  This  is  done  in  appendix  D. 
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The  characteristic  function  method  presented  in  appendix  A  provided  a 
powerful  tool  to  determine  the  moments  for  the  many  processes  included 
in  the  dissertation. 
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II  MODELING  OF  PHASE  AND  FREQUENCY  INSTABILITY 


2.0.  Introduction 

As  indicated  earlier,  by  frequency  instability  of  a  precision 
frequency  source  we  mean  unwanted  frequency  departures  from  a  nominal 
value.  Frequency  stability  describes  the  degree  to  which  a  source 
produces  a  constant  frequency  over  a  specified  time  interval.  In  this 
chapter,  we  consider  the  characterization  of  this  time-dependent 
phenomenon. 

A  certain  class  of  non-stationary  as  well  as  stationary  phase 
random  processes  will  be  identified  through  their  probability  density 
functions  and  their  moments,  particularly  their  covariance  matrices. 
This  will  allow  determination  of  the  auto-correlation  function  of  the  rf 
signal,  i.e.  the  oscillator  output  sinusoidal  signal.  It  will  be 
shown,  for  such  a  class  of  phase  random  process,  that  the  rf  signal  is 
stationary  In  auto-correlation  function  despite  the  basically 
non-stationary  phase  driving  force.  In  chapter  three,  we  will  establish 
conditions  under  which  the  rf  signal  is  wide  sinse  stationary. 

The  method  used  here  provides  a  different  way  for  determining  the 
time  and  frequency  behavior  of  oscillators.  Most  prior  work  has  dealt 
with  the  non-stationary  aspects  of  phase  through  the  use  of  finite  time 
averages  (where  time  shifted  auto-correlation  functions  were  not 
addressed).  By  commencing  with  the  phase  covariance  properties,  the 
direct  use  of  non-stationary  phase  instability  leads  nonetheless  to 
stationarity  of  the  rf  signal  auto-correlation  function.  In  a  later 
chapter,  the  spectral  spreading  (dispersion)  due  to  phase  instability 
will  easily  follow  from  the  material  presented  in  this  chapter. 
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2.1.  Oscillator  RF  Signal 


The  low  pass  equivalent,  b(u),  of  the  ideal  oscillator  output 
waveform,  say  r(u),  is  defined  through 

r(  u  )  -  Re  {  b(  u  )  exp(  jw  u  )  } 

o 

where  w  is  the  oscillator  frequency.  b(u)  may  be  written 
o 

j*(u) 

b(  u  )  »  B  e 

where  B  is  some  arbitrary  amplitude  and  0(u)  is  a  deterministic  phase 
function.  For  a  monochromatic  oscillator  0(u)  «*  a  constant.  It  is 
this  low  pass  equivalent,  b(u),  which  we  will  consider  throughout  this 
dissertation  and  will  refer  to  it  as  the  "oscillator  rf  signal”  in 
keeping  consistant  with  the  literature  [2.1],  [2.2].  For  non-ideal 
oscillators  found  in  practice,  there  will  be  unwanted  departures  from 
the  nominal  amplitude  and  frequency.  We  may  express  the  output 
oscillator  waveform  as 

j[0  +  z(u  )] 

b(u  )  *  [B  +  c(u  )]  e 

where  c(u)  denotes  the  amplitude  noise  and  z(u)  denotes  the  phase  noise. 
Both  of  these  are  random  processes.  Since  most  oscillators  are 
amplitude  stabilized  [2.3],  we  shall  neglect  the  amplitude  noise  and, 
without  loss  of  generality,  assume  B  to  be  unity.  Also,  since  |/  is  a 
deterministic  phase,  we  shall  consider  the  following  simplified 
oscillator  signal  model. 


2 


s 


j  z(u) 


(2.1) 


b(  u  )  “  e 


In  the  above  oscillator  rf  signal,  the  phase  random  process,  z(u), 
generates  the  random  phase  fluctuations.  We  shall  focus  our  attention 
on  the  modeling  of  this  phase  error  process,  z(u).  The  relationship 
between  the  frequency  noise  and  the  associated  phase  random  process  has 
been  discussed  in  [2.1]— [2.4] .  Basically,  a  phase  random  process  is  the 
integral  of  the  frequency  random  process.  Translating  variances  and 
spectrum  from  frequency  to  the  phase  random  process  must  carefully 
account  for  this  integral  relationship. 

Traditionally,  the  phase  error  process  has  been  assumed  to  be  a 
stationary  process  [2. 4], [2. 5].  The  reason  has  been  the  simplicity  of 
such  an  assumption.  However  this  assumption  has  limitations  as  pointed 
out  in  [2.4].  Certain  idealized  models,  [2.3],  for  phase  and  frequency 
random  processes  impose  serious  limitations  on  their  application  to 
stability  analysis  by  excluding  a  whole  class  of  phase  random  processes 
from  such  consideration.  Examples  include  the  random  walk  phase  and 
random  walk  frequency  processes. 

The  Wiener  (or  BM:  Brownian  Motion,  random  walk)  phase  random 
process  has  no  derivative  and  is  furthermore  non-stationary ;  in  other 
words  it  has  no  instantaneous  frequency  and  cannot  be  represented  by  a 
power  spectral  density.  In  this  dissertation,  these  difficulties  are 
overcome  by  assuming  a  finite  bandwidth  oscillator,  i.e.  a  "near  white” 
decorrelation  time  and  by  focussing  attention  on  the  output  signal, 
b(u  ). 


Since  we  are  dealing  with  the  time  and  frequency  domain  properties 
of  the  oscillator  rf  signal  itself,  we  find  it  appropriate  to  focus 
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attention  on  the  sinusoidal  output  signal,  b(u),  and  its  power  spectral 
density,  S(w),  [2.4].  Knowledge  of  the  power  spectral  density  is  of 
prime  importance  in  applications  such  as  radar  Doppler  processing, 
communications,  frequency  synthesis,  and  spectroscopy.  When  the  AH 
noise  is  negligible  and  the  phase  instability  is  much  smaller  than  one 
radian,  then  the  low  modulation  index  approximation  can  be  used,  i.e. 
the  spectral  density  of  the  waveform,  b(u),  is  approximated  by  the  noise 
sidebands  of  the  phase,  z(u).  However,  the  complete  formula  is  often 
needed  to  forcast  spectral  purity  when  very  high  ratio  frequency 
multiplication  is  implemented,  [2. 4], [2. 6]. 

We  also  make  use  of  the  fact  that  for  a  class  of  non-statlonary 
phase  random  processes,  z(u),  the  sinusoidal  transformation  (or 
operation)  yields  a  random  process  which  is  stationary  in 
auto-correlation  function.  Specifically  we  obtain  stationary 
auto-correlation  properties  for  the  signal,  b(u),  for  the  case  when 
first  and  higher  order  increments  of  z(u)  are  stationary.  In 
chapter  III,  we  determine  conditions  for  a  general  class  of 
non-s tat ionary  phase  random  processes  which  lead  to  wide  sense 
stationary  rf  signals  through  the  above  sinusoidal  operation.  In  a 
later  chapter,  we  will  find  that  the  parameters  within  the  argument  of 
its  auto-correlation  function  have  a  one-to-one  correspondence  with  the 
power  spectral  density.  We  will  make  use  of  this  condition. 
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2*2.  Modeling  of  Phase  and  Frequency  Components 


Before  we  proceed  with  the  modeling  we  shall  describe  the  properties 
l 

of  the  Wiener  process  .  From  Ross  [2.7]  the  Wiener  process,  say  x(s), 
has  the  following  properties: 

A  stochastic  process  [x(s),  s  0]  is  said  to  be  a  Wiener  process 

x 

with  mean  drift  coefficient,  jt,  and  variance  drift  coefficient,  <r  ,  if; 

i.  x( 0)  -  0 

ii.  (x(s),  8^0}  has  stationary  independent  Increments 

iii.  For  every  s  >  0,  x( s )  is  normally  distributed  with  mean, 
(aXs),  and  variance,  s. 


Since  we  want  the  time  when  the  observation  is  initiated  to  be 
arbitrary,  our  model  will  be  modified  as  follows: 


We  define  a  random  walk  phase  process  to  contain  a  phase  random 

variable,  x(u)  -  X,  corresponding  to  the  condition  at  the  time, 

u(~o©  <  u  <cw),  when  the  observation  is  initiated.  This  random 

variable,  X,  is  assumed  to  have  a  Gaussian  probability  density  function 

where 


E{  X  } 

and 


’  ar  {  X  } 


an  arbitrary  constant 
2 

c  <f  -an  arbitrary  constant. 
1  1 


l 

The  Wiener  process  is  sometimes  refered  to  Brownian  motion 
or  the  Gaussian  random  walk  process.  We  will  frequently  use  these 
designations  throughout  this  dissertation. 


We  add  to  this  random  variable,  an  independent  Wiener  process,  say  x  (  t ), 

o 

where  x  (0)  A  0,  the  zero  mean  drift,  ^  ■  0  and  the  drift  variance  is 
2  o 

<f  •  Furthermore,  our  process  is  defined  to  run  for  negative  time  as 
1 

well  as  positive  time,  tracing  mirror  images  of  paths  with  equal 
probability  from  time,  u.  Thus,  t  will  be  replaced  by  |t|.  From,  the 
independence  assumption, 

x(  u+t  )  =  X  +  x  (  t  ) 
o 

2  2 

implies  var{x(u+t)}  ■  var{  X  }  +  var{x  ( t )}  ■  c  <T  +  <J'|t| 

o  111 

We  shall  use  the  term,  x(u),  to  denote  the  random  walk  phase  process 
throughout  this  dissertation.  Next  we  briefly  describe  the  various 
types  of  phase  instability  models  appearing  in  the  literature. 


Generally,  phase  random  processes,  z(u)  (not  b(u)),  are  modeled  as 
power  law  spectral  densities.  They  have  been  described  in  the 
literature  and  have  won  wide  acceptance  [2.4].  The  phase  instabilities 
are  modeled  according  to  the  rate  at  which  the  power  spectral  density  of 
the  phase,  z(u)  (or  frequency,  z(u))  random  process  decays  with  the 
transformed  variable,  f.  Hence,  if  S^( f )  is  the  power  spectral  density 
of  the  frequency  random  process,  z(u)  (the  derivative  of  the  phase 


°  <  f  <  f 

h 


f  >  f 

h 

where  f  is  an  upper  cutoff  frequency, 
h 

The  models  for  various  values  of  a  ■  -2,-1,  0,  1,  2  correspond  to  random 


random  process,  z(u)),  then 

+2 


S.(f  ) 

z 
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walk  frequency,  flicker  frequency,  white  frequency/random  walk  phase, 
flicker  phase,  and  white  phase  respectively.  The  flicker  models,  i.e. 
a  ■  -1,  1,  will  not  be  treated  here. 

The  models  we  will  use  tend  to  have  certain  physical  origins, 
[2.4],  which  are  as  follows: 

a)  White  phase  noise,  y(u),  is  usually  due  to  white  noise  sources 
external  to  the  oscillator  loop. 

b)  Random  walk  phase  (white  frequency)  noise,  x(u),  usually 
arises  from  additive  white  noise  internal  to  the  oscillator 
loop. 

c)  The  phase  random  process,  v(u),  resulting  from  random  walk 
frequency  noise,  v(u),  is  usually  related  to  the  oscillator 
environment  such  as  temperature,  vibration  and  shocks. 


Specifically  we  assume  that  the  phase  process  is  composed  of 
non' stationary  components,  v(u)  +  x(u)  and  stationary  components,  y(u), 
so  that  the  phase,  z(u),  is  modeled  as 

z(u)«v(u)  +  x(u)  +  y(u)  (2.2) 

• 

where  v(u)and  x(u)  are  zero  drift  mean  Wiener  processes,  characterized 
as  having  stationary  independent  increments  and  which  gives  rise  to  the 
non-s tat ionary  random  walk  frequency  and  phase  [2. 4], [2. 5]  instability 
respectively.  y(u)  is  a  stationary  Gaussian  process  which  gives  rise  to 
the  white  x  [2. 4], [2. 5],  instability.  Since  we  are  considering  finite 
bandwidth  oscillators,  the  random  process,  y(u),  decorrelates  over  some 
interval  -To/2  <  t  <  To/2.  We  will  refer  to  this  process  as  "near 
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white".  (Note  notational  change  from  {2. 4],  [2. 5];  in  our  work  v(u), 
x(u)  and  y(u)  are  each  phase  functions). 

As  indicated  earlier,  white  phase  noise,  y(u),  and  random  walk  phase 
noise,  x(u),  are  caused  by  additive  white  noise  both  external  and  internal 
to  the  oscillator  loop  where  as  the  random  walk  frequency,  v(u),  is  due 
to  the  oscillator  environment.  In  this  dissertation,  we  shall  develop 
models  for  the  internal/external  white  noise  induced  instability 
separately  from  the  instability  induced  by  the  oscillator  environment. 


2.2.1.  Model  for  White  Phase  and  Random  Walk  Phase 

First  we  consider  the  models  where  the  phase  random  process  is  the 
superposition  of  the  stationary  white  phase  process,  y(u),  and  the 
random  walk  phase  process,  x(u).  The  means  of  these  phase  random 
processes  can  be  arbitrary.  However,  without  loss  of  generality,  we 

1^ 

assume  that  the  means  are  zero  . 

Since  oscillators  have  a  finite  bandwidth,  let  us  assume  a 
triangular  rather  than  an  impulse  auto-correlation  function  for  the 
stationary  white  phase  process  y(u)  as  shown  in  figure  2.1a.  This  is 
consistent  with  the  assumption  that  y(u)  gives  rise  to  the  "near  white" 
instability  with  the  area  under  the  triangular  (impulse  like)  part  of 
the  auto-correlation  function  being  equal  to 

_ 

When  evaluating  the  auto-correlation  function  of  b(u)  it  will 
become  obvious  that  mean  values  of  x(u),  y(u)  and  v(u),  c.f.  equation 
(2.2),  cancel  out  of  the  exponential  expression.  Therefore,  the  mean 
values  of  the  phase  random  processes  do  not  affect  the  auto-correlation 
function  of  b(u). 
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TRIANGULAR  AUTO-CORRELATION  FUNCTION 
STATIONARY  NEAR  WHITE  PHASE  STABILITY 


RANDOM  WALK  COMPONENT  CORRELATION 

OF  STABILITY  FACTOR  X(u) 

(b) 


FIGURE  2.1.  STATIONARY  AND  NON-STATIONARY  STABILITY 


2 

<r  . 

o 


As  seen  in  figure  2.1a,  we  have: 


2 

Var{y(u)}  =  Var{y(u+t)}  ■  c  o' 

o  o 


(2.3a) 


Cov{y(u+t)  y(u)}  -  E{y(u+t)  y(u)} 


(2.3b) 


2  2  To 

Cov{y(u+t  )  y(u)}  «  <f  (c - |t|)  :  |t|  <  — 

o  o  To  2 

2  To 

-  <T  (c  -  1)  :  |  t|  >  —  (2.3c) 

o  o  2 


l 

where  c  is  arbitrary  . 


I 

The  phase  random  processes,  x(u),  y(u)  and  v(u)  are  functions  of 
time,  u.  However,  it  will  become  obvious  that  in  the  case  of  the  rf 
signal,  b(u),  its  auto-correlation  function  depends  only  on  the  phase 
change  from  time,  u,  to  time,  u+t,  and  not  on  the  value  or  statistics 
at  time,  u,  -ot>  <  u  <  a©.  We  will  therefore  associate  arbitrary  values, 
c  ,  c  and  c  to  the  phase  statistics  at  time,  u. 
o  1  3 
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The  variance  of  the  random  walk  component  is  shown  in  figure  2.1b. 
From  the  definition  of  our  random  walk  process,  with  u  being  the 


reference  time,  the  following  can  be  written  [2.8]  -  [2.11] . 

2 

Var{x(u)}  *  c  cr  (2.4a) 

1  1 

2 

Var  {x(  u+t  )}  ■  (  c  +1  1 1  )  <r  (2.4b) 

1  1 

2 

Cov [x( u+t  )  x( u  )}  =  c  c T  (2.4c) 

1  1 

2  I 

where  <T  is  the  drift  variance  of  the  random  walk  (Wiener)  component 

1 

and  c  is  arbitrary.  This  assumes  that  the  variance  of  x(u)  at  some 
1 

reference  point,  u,  is  arbitrary  and  that  stationary  zero  mean 
independent  Gaussian  increments  (whose  variance  is  proportional  to  t ) 
are  additive  to  the  process  with  increasing  or  decreasing  t.  The 
Independent  increment  property  can  be  used  to  derive  equation  (2.4c)  as 
follows  [2.7]; 

Cov {x( u+t)  x(u)}  ■  E{[x(u+t)  -  x(u)  +  x(u)]  [x(u)]i 

*  E[[x(u+t)  -  x(u)]  [x(u)]}  +  E{[x(u)j  [x(u)]} 

2  2 
■  0  +  c  <f“c  (T 
11  11 


which  is  just  the  variance  of  the  random  variable,  x(u),  defined  at  the 

beginning  of  this  section  at  time  =  u,  i.e.  t  =  0.  Later  we  shall  see 

f  2  2 

The  variances  in  this  paper  denoted  by  o'  and  c r  correspond 
to  h  and  h  respectively  in  [2.4]  ,  [2.5]  .  o  1 

2  o 


2  2 

that  the  arbitrary  variance  c  tf"  and  c  o'  which  we  assign  to  x(u)  and 

1  1  o  o 

y(u)  respectively  will  not  affect  the  auto-correlation  function  of  the 

oscillator  rf  signal,  b(u).  When  determining  the  auto-correlation  and 

spectral  properties  for  the  case  when  the  envelope  waveform,  f  (u), 

o 

is  constant,  i.e.  a  cw,  the  oscillator  output  itself,  we  use  the  complete 

description  of  equations  (2.3)  and  (2.4).  For  the  case  when  f  (u)  is  a 

o 

train  of  pulses  we  can  assume  that  the  phase,  z(u),  is  constant  within 
each  of  the  small  intervals  of  d(u-iT)  denoted  by  /\u  .  In  doing  this, 
we  are  assuming  that  the  pulse  width,  Au>  is  much  smaller  than  the 
interpulse  period,  T,  i.e.  /\u  <<  T. 

There  are  two  major  approaches  for  the  frequency  stability 
characterization.  One  deals  with  the  characterization  in  the  time 
domain  and  the  other  with  the  characterization  in  the  frequency  domain. 
Translations  between  the  two  domains  have  also  been  described  in  the 
literature  [2. 4], [2. 5].  For  the  most  part,  this  has  been  limited  to 
time  averages  of  the  phase /frequency  random  processes  including 
different  ways  of  determining  their  variances.  Equivalent  spectral 
domain  representations  of  the  phase /frequency  random  processes  include 
the  use  of  transformed  window  functions  acting  on  the  phase  process 
itself  resulting  in  sine-x-over-x  convolutions  in  the  spectral  domain. 
In  this  chapter,  we  shall  focus  on  the  details  of  the  time  domain  rf 
signal  auto-correlation  function.  This  will  set  the  stage  for  spectral 
domain  characterization.  The  relationship  between  the  two  domains  will 
thus  be  established  in  the  subsequent  chapters.  The  power  spectral 
density  of  the  complete  oscillator  signal  output,  b(u),  as  affected  by 
the  phase  random  process,  z(u),  will  be  determined.  Characterization  in 
the  frequency  domain  is  important  both  for  theoretical  reasons  and  for 
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application  purposes.  Spectral  spreading  is  a  primary  specification  for 
engineers  in  spectroscopy  and  radar  and  it  is  therefore  a  widely 
accepted  chatacterization  of  oscillator  stability. 

The  oscillator  waveform  was  written  in  equation  (2.1)  as  a  low  pass 
equivalent  to  a  sinusoidal  random  process  containing  the  phase  process, 
z(u),  as  follows 

j  z(u) 

b(u)  =  e  (2.1) 


The  modulating  envelope  waveform,  f  (u),  may  be  written 

o 

f  (u)  ±  c 
o 


for  the  cw  case  and 


N/2-1 


I  T  v— 

f  (u )  A  c| - \  d(  u 

o  I  Au  , 

\j  i— N/2 


I  i  +  —  ]  T  ) 
2 


(2.5a) 


(2.5b) 


for  the  pulse  train  of  N  (even)  pulses  where  the  interval,  Au»  an^ 
the  function,  d(u),  are  defined  through 


From  equations  (1.2),  (1.3)  and  (1.4)  we  may  evaluate  the  average 

auto-correlation  function,  R  ( t ),  of  the  system  output  waveform,  f(u), 

f 

by  forming  the  product  of  the  auto— correlation  function  of  the 


oscillator  output  waveform,  equation  (2.1),  and  the  time  average 


auto-correlation  function  of  the  modulating  envelope  waveform,  equation 
(2.5).  The  auto-correlation  function  of  the  oscillator  output  waveform, 

equation  (2.1),  is  [2.12] , [2.13] , [2.14] 


V 


R  (t)  ^  E{b(u  +  t)  b  (u)} 
b 


j[z(u+t)]  -j[z(u)] 

=  E{e  e  } 


j[z(u+t)  -  z(u ) ] 

-  E{e  } 


j[  x(u+t)  +  y(u+t)  -  x (u )  -  y(u)] 

-  E{e  } 


The  argument  of  the  exponential  can  be  expressed  in  vector  notation 
by  defining 


and 


x(u+t) 

y(u+t) 

+ 

X(u) 

y(u) 

j  [  P  , 
1 


P  I 
2 


(2.6) 


In  Appendix  A,  we  illustrate  the  use  of  characteristic  functions  for  the 

evaluation  of  higher  order  moments  of  the  rf  signal,  b(u).  Thus  the 

auto-correlation  function  of  the  oscillator  rf  signal  is  written  using 

the  procedure  outlined  in  Appendix  A.  We  have 

*  jz(u+t)-jz(u) 

R(t)^E(  b(u+t)b  (u)}  -  E  (  e  ) 

b 

which  in  vector  form  is 
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z_  2_  1 

R  (t)  -  E{  e  }  *  F  (p)  -  exp{ - (  m  +  m  -  2m  )}  (2.7) 

b  z  2  11  22  12 


114 


which  is  the  characteristic  function,  F,  of  z  evaluated  at 


fit 


L  " 


L-U 


and  where  m  is  the  ij  th  element  of  the  covariance  matrix,  M 

ij  1 

In  vector  notation  for  z(u),  the  joint  density  function  of  a 
Gaussian  random  process,  z_t  is  [2.7],  [2.14] 


P  {z} 


2  IT  |M 


1  T  -1 

- exp  {-  —  z  M  z } 

1/2  2  z 


(2.8) 


where  the  superscript,  T,  denotes  transpose.  The  above  covariance 

matrix,  M  ,  which  is  derived  from  equation  (2.3)  and  equation  (2.4) 
z_ 

through  Appendix  B  is  given  by 


M  -  M  +  M 
—  *  i 


r  2  2 

I  ff(c  +  lt|)  +  <T  c 


2  2  2 

<r  c  +  cr  (c  -  —I  tl  ) 
11  o  o  To 


2  2  2 

<T  c  +  <r  (c - 1  tl  ) 

11  o  o  To- 


2  2 

<r  c  +  <T  c 
11  o  o 


(2.9a) 


for  Itl  <— -  and 
2 
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r  2  2 

I  (T  (c  +  I  t|  )  +  cr  c 
111  o  c 


2  2 

o-  c  +  <r  (c  -  1) 
11  0  o 


2  2  1 
<J*  c  +  <J*(c-l)l 
1  1  o  o  I 


2  2 

(T  C  +  <J*  C 
11  o  o 


(2.9b) 


To 

for  | t I  >  — ~ 

2 

We  now  find  the  auto-correlation  function  of  the  oscillator  rf 

signal,  R  (  t  ),  using  equation  (2.7)  and  the  value  of  £  given  above, 
b 

We  obtain: 

T 

j[z(u+t  )-z(u  )]  ££ 

R  (t)  =  E{e  }  -  E{e  } 

b 


exp  {  £  £  }  p  {£}  d£ 
z 


which  gives 


1  2 


R  (  t  )  =  exp  {-  — ( <r  1 1 1  +  2  <T - I  1 1  )}  :  I  1 1  <  — 


2  1 


o  To 


1  2  2 

exp{-  —  (cr  |  tl  +  2  <r  )} 
2  1  o 


:  I  tl  >  --  (2.10) 

2 


The  logarithm  of  R  (  t  )  is  shown  in  figure  2.2. 
b 
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FIGURE  2.2  WHITE  PHASE  AND  RANDOM  WALK  PHASE,  LOGARITHM  OF 

AUTO- CORRELATION,  R  (  t )  OF  OSCILLATOR  WAVEFORM  b(u) 

b 


It  Is  well  known  [2.14],  [2.15]  that  the  Fourier  Transform  of  the 
auto-correlation  function  gives  the  power  spectral  density  when  the 
process  is  a  stationary  random  process.  We  establish  conditions  for 
stationarity  and  ergodlcity  in  chapter  III  and  proceed  to  evaluate  power 
spectral  densities  as  well  as  average  power/energy  spectrums  for 
different  modulating  envelopes  in  chapter  IV. 
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2.2.2.  Model  for  Random  Walk  Frequency 


As  mentioned  previously  random  walk  frequency  is  yet  another 
instability  factor  identified  in  oscillators.  This  is  caused  by  the 
environment  in  which  the  oscillator  must  operate.  The  random  walk 
frequency  model  is  characterized  by 


2 

var  {v(  u  )}  =  c  O' 
2  2 


(2.11a) 


2 

var  {v(  u+t  )}-(c  +  |t|)<r  (2.11b) 

2  2 


2 

cov{v(u+t)  v(  u  )}  =  c  &  (2.11c) 

2  2 

2  _1 
where  c  cr  is  the  variance  of  the  frequency  offset  due  to  the  effects 
2  2 

of  long  term  drift  over  time,  u.  The  frequency  variance  term  due  to  the 
drift  variance  coefficient  has  a  linear  slope  as  shown  in  figure  2.3a. 
The  variance  of  the  resulting  phase  has  a  third  power  relationship  with 
time  as  shown  in  figure  2.3b  and  proved  in  Appendix  C.  The  properties 
of  the  rf  signal  associated  with  this  process  can  be  investigated  in 
much  the  same  manner  as  in  the  case  of  the  above  random  walk  phase 
(white  frequency)  component.  We  will  first  determine  the  covariance 
matrix  of  the  phase  random  process  associated  with  the  random  walk 
frequency. 


I 

Unlike  the  arbitrary  constants,  c  ,  c  and  c  (associated  with 

o  1  3 

the  phase  statistics  for  any  u;  -ar>  <  u  <  oo  )  which  have  no  effect  on 
the  auto-correlation  function  of  the  rf  signal,  the  constant,  c  , 

2 

(associated  with  the  long  term  frequency  linear  drift  over  all 
u;  -at>  <u  <oo)  does  affect  the  auto-correlation  function  over  the 
short  term  time  interval,  denoted  by  t.  It  will  be  shown  later  that 
this  term  can  be  treated  as  a  constant  over  time,  t,  if  the  change 
with  respect  to  time,  u,  is  very  slow. 


-  36  - 


Using  this  covariance  matrix,  the  auto-corrclation  function  of  the 
oscillator  signal  containing  the  random  walk  frequency  factor  will  be 
derived  next  through  the  characteristic  function  approach  given  in 
Appendix  A.  Third,  the  power  spectral  density  is  obtained  in  chapter  IV 
by  taking  the  Fourier  transform  of  the  auto-correlation  function. 
Spectral  spreading  due  to  the  random  walk  frequency  component  can  then 
be  determined. 

In  much  the  same  fashion  as  in  Section  2.2.1,  we  find  the 

auto-covariance  matrix  of  the  process,  v(u),  which  results  from  the 

random  walk  frequency,  v(u).  The  phase  random  process  may  consist  of  an 

arbitrary  phase  random  variable,  v(u  ),  including  accumulated  drift 

o 

effects  [2.1], [2. 2]  up  to  some  time,  u  ,  plus  the  integral  of  the 

.  o 

frequency  random  process,  v(u). 

u  +u 
°(  • 

v(u )  *  v(u  )  +  \  v(h )  dh  (2.12) 

u 

o 


This  is  a  non-s tat ionary  random  process  for  which  the  auto-covariance 
matrix  can  be  derived  as 


From  (2.13)  we  see  that  we  must  find 


-  38  - 


* 

var{v(u)}  «  E{  v(u)  v  (u)} 


* 

var{v(u+t)}  *  E{v(u+t)  v  (u+t)} 


cov  {v(  u+t  )v  (  u  )}  ■  E  {v(  u+t  )  v  (  u  )} 

The  above  expectations  have  been  evaluated  in  Appendix  C  whereby  the 
covariance  matrix  for  the  phase  random  processes,  v(u),  is  given  by 


2 


3 

2  Itl 

c  +c  |t|  + - 

3  2  3 


<r 

2 


c 

3 


(2.14) 


2  2 
where  C  is  the  random  walk  frequency  drift  variance,  c  <r  is  the 
2  2  2 

2 

variance  of  the  frequency  due  to  long  term  drift,  and  c  <r  is  arbitrary. 

3  2 


Since  the  process,  v(u),  is  Gaussian,  we  can  use  Appendix  A, 

equation  (A. 3b),  to  evaluate  the  auto-correlation  function,  R  ( t ),  of  an 

b 

oscillator  RF  signal,  b(u).  This  is  a  random  process,  whose 

auto-correlation  function  is  invariant  with  time,  u,  and  dependent  only 

on  the  shift  variable,  t.  This  property  of  R  ( t ),  based  on  the 

b 

covariance  matrix  of  equation  (2.14),  will  be  established  In  greater 
detail  through  Theorem  1  and  Theorem  II  in  chapter  III. 


R  (t) 
b 


r  2  2  3 

I  <T  (  3c  Itl  +  Itl  ) 


(2.15) 
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2.2.3.  Phase  and  Frequency  Linear  Drift 


? 


Frequency  drift  due  to  aging  of  the  oscillator  results  in  an 
extremely  slow  variation  of  the  frequency  process  with  respect  to  time. 
This  can  be  approximated  by  a  linear  term  about  some  local  point  [2.1] , 
[2.2] .  Let  us  define  the  total  phase  random  process  as  follows: 

•  • 

z(u)*a  +  bu  +  z(u)  (2.16) 

T 

We  can  proceed  to  find  the  auto-correlation  function  of  the  rf  signal. 
Assuming  a  and  b  are  random  variables  and  z(u)  is  an  independent,  zero 
mean,  random  process,  all  Gaussian,  where  z(u)  contains  the  processes 
as  described  previously  in  equation  (2.2),  we  determine  the  phase  random 
process  of  the  entire  signal. 

b  2 

z(u)*«c+au  +  -u+z(u)  (2.17) 

T  2 


The  covariance  matrix  is  defined  by 


T 

M  ^  E(  z  z 
z^  T  T 

T 


2  2  2  2  2  2 
Now  <f  *E(c  ),  <f  ^E(a  ),  cr  ^  E(  b  ). 
cab 


Also,  time,  u  and  t,  are  non-random,  i.e. 


*  A  *  t  . 

E(u  t  )  *  u  t:  k,/  *  0,  1,  2,  .... 
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Furthermore,  the  covariance  matrix  of  z_  has  been  given  for  each  of  its 

independent  components,  x,  jr,  and  z_,  through  equation  (2.9)  and  equation 

(2.14).  Then  the  covariance  matrix  of  z_  ,  above  is  easily  written: 

T 


f  2  2 

I  <r  c 

{2  2  2  b  4  22  b22 

I  O'  4tT  ( u+t  )  +  —  (  u+t  )  4m  o'  +<f  u(  u4-t  )  4-  —  u  (  u4-t  )  4m 

lea  4  11  c  a  4  12 


.  2 
<T 

2  2  b  2  2 

<f  4tr*  u(  u4-t  )  4-  —  u  (  u4-t  )  4m 
c  a  4  21 


2 

O' 

2  2  2  b  4 

O'  Hf  u  4-  —  u  4m 
c  a  4  22 


(2.18) 

where  m  la  the  lj  th  element  of  the  covariance  matrix,  M  ,  defined 

ij  L 

in  Appendix  A,  and  c  is  an  independent  Gaussian  random  variable.  The 
auto-correlation  function  of  the  rf  signal  is  found  by  taking  the 
characteristic  function 


1  T 

R(u,t)«F  (jO  ■  exp  {-  -  £  M  2.  { 
b  z_  2  z_ 

T  T 


1  2  2  2  2  2 

-  exp{ - (  o'  t  4-  o'  t  [ 2u4-t ]  4-  J  m  )}  (  2.19 ) 

2  a  b  ij 


T 

where  £  ■  [  41,  -1  J  and  the  appropriate  m  are  obtained  from 

ij 

equation  (2.18)  through  the  formula  in  equation  (2.7). 

2 

The  variance  of  the  phase  accumulated  up  to  time  u  is  o'  and  the 

c 

2  2  2 

variance  of  the  frequency  accumulated  up  to  time  u  is  o'  4-  o'  u  . 

a  b 


A  further  consideration  of  equation  (2.19)  suggests  that  the  first  and 
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second  order  drift  terms  (first  and  second  terms  in  the  exponent  of 
equation  (2.19))  can  be  approximated  as  follows: 


k 


y 


r* 


m 


i 


22  22  2  222  2 
(<r  +  <r  [2u+t]  )t  2(<r+f2u0']  l<r(u,t)t  (2.20) 
a  b  a  b  drift 


for  0  <  t  <<  u.  The  drift  variance  resulting  from  this  approximation  is 
2 

lumped  into  c  o'  for  the  offset  frequency  in  the  treatment  of  the  random 
2  2 

walk  frequency  instability  model  in  Section  2.2.2.  In  general,  it  is  a 
function  of  u  and  t  where  u  gives  rise  to  long  term  drift  and  t  is  the 
short  term  shift  variable  which  has  a  negligible  effect  on  the  drift 
variance.  From  equation  (2.16),  it  has  an  assumed  Gaussian 
distribution. 


It  should  be  noted  that  the  expressions  in  the  exponent  of  the 
auto-correlation  functions,  equations  (2.10),  (2.15),  and  (2.19) 

represent  the  first  increment  structure  function  with  tl  *=t2  “  t. 
Structure  functions  used  by  Lindsey  and  Lewis  [2.1J  were  first 
introduced  by  Kolmogorov  [2.16].  The  first  increment  structure  function 
of  a  phase  random  process  (instability)  Is  defined  by 


D  ( u  ;t  ,t  )^E  {  [  z(t+u)-z(u)]  [  z(t  +  u)-z(u)j  } 
z  1  2  1  2 


The  structure  functions  are  useful  in  obtaining  phase  stability  measures 
[2.1]  as  we  shall  demonstrate  in  the  following  section. 


2.3.  Frequency  Stability  Measures 

Thus  far  in  this  chapter,  we  have  considered  a  non-stationary  model 
for  the  phase  random  process  from  which  the  covariance  matrices  and 
auto-correlation  functions  of  the  rf  signal  followed.  Now  we  consider 
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phase  stability  measures  to  describe  precision  frequency  sources  and  how 


their  frequency  instability  affects  system  performance. 


An  appropriate  measure  to  use  in  conjunction  with  phase  instability 
models  discussed  in  this  chapter  is  based  on  the  work,  of  Lindsey  and 
Lewis  [2.1]  and  Lindsey,  et  al  [2.2]  .  This  involves  the  elements  of  the 
phase  covariance  matrix  developed  throughout  this  dissertation.  The 
measures  are  dimensioned  in  phase  standard  deviation  and  they  provide  a 
quantitative  description  of  the  phase  deviation  with  respect  to  time. 
For  all  the  models  discussed  in  this  chapter,  the  covariance  matrices 
have  been  derived  and  the  frequency  stability  measure  can  be  computed  in 
a  straight  forward  manner  through  use  of  the  structure  function.  From 
[2.1] ,  [2.2]  we  have 


dz(t  ) 


i  D  (t) 
\1  * 


(2.21a) 


Generally,  in  terms  of  the  covariance  matrix  elements,  m  ,  of  M  , 

ij  z 

given  in  equations  (2.9),  (2.14)  and  (2.18),  we  have 


dz(  t  ) 


m 

N 


+m  -2m 
11  22  12 


i*“ - -  - 

-2  In  R  (  t  )  =  -2  In 

\|  b  \1 


F  (l,  -1) 

z 


(2.21b) 


This  is  the  square  root  of  the  exponent  in  the  auto-correlation 
function;  given  that  the  reference  time  for  all  processes  is  u.  In 
Table  2.1,  we  provide  the  square  of  the  frequency  stability  measures  of 
equation  (2.21)  for  each  of  the  processes  considered  in  this  chapter. 


Next,  we  investigate  stochastic  properties  of  the  rf  signal,  b(u), 
as  affected  by  our  phase  instability  models.  Stationarity  and 
ergodicity  in  the  mean,  auto-correlation  function  and  power  spectral 
density  are  all  considered. 
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TABLE  2.1  MEASURES  OF  PHASE  INSTABILITY 


PROCESS 

-2  in  F  in) r  D  1 1) 

2  —  Z 

WHITE  PHASE 

fS'ti 

STARTING  PHASE 

0 

RANDOM  WALK  PHASE 

«•,*  Itl 

STARTING  FREQUENCY 

o’  a  Itl2 

RANDOM  WALK  FREQUENCY 

i  z  i.i* 

3*2  Ul 

LINEAR  DRIFT  FREQUENCY 

(2u<rb)2|t|2 

UA  LONG  TERM  P  s  [  j."  ij 

t  A  SHORT  TERM 
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III.  STATIONARITY  AND  ERGODICITY  OF  THE  OSCILLATOR  RF  SIGNAL 


3.0 .  Introduction 

In  the  previous  chapter,  we  discussed  three  models  for  the  phase 
random  process,  i.e.  "near  white"  phase,  random  walk  phase,  and  random 
walk  frequency  random  processes.  In  this  chapter,  we  will  study  the 
stationarity  and  ergodicity  properties  [3.1],  [3.2]  of  the  rf  signal 
driven  by  the  above  three  phase  random  processes.  Conditions  will  be 
established  for  the  stationarity  and  ergodicity  in  the  mean, 
auto-correlation  function,  and  power  spectral  density.  These  results 
will  be  useful  in  the  justification  of  what  is  to  follow  in  chapter 
four,  i.e.,  the  determination  of  the  power  spectral  density  through  the 
Fourier  transform  of  the  auto-correlation  functions.  First,  we  consider 
the  stationarity  properties. 

3.1.  Stationarity 

We  now  examine  conditions  for  stationarity  of  the  rf  signal. 
Rutman  [3.3],  and  Boileau  and  Picinbono  [3.4]  and  others  have  called 
attention  to  the  fact  that  one  should  not  misuse  the  stationarity 
assumption  or  the  correlation  function,  R(  t  ),  that  has  been  widely  used 
in  the  frequency  stability  literature.  Although  it  is  convenient  to 
assume  stationarity,  one  must  check  that  this  is  not  in  conflict  either 
with  other  parts  of  the  model  or  with  some  physical  arguments  [3.3].  As 
already  pointed  out,  our  models  for  the  phase  random  process,  z(u),  are 
generally  non-stationary .  Since  our  attention  is  being  focused 
primarily  on  the  oscillator  rf  signal. 
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j  z(u  ) 

b(u)^e  (3.1) 


we  will  proceed  to  establish  the  stationarity  of  b(u). 


Definition  1: 


A  random  process  which  satisfies 


E{b(u  )}  =  \  (  3.2a) 

* 

and  R(t)^E{b(  u+t )  b  ( u )  }  (3.2b) 

b 


for  all  u  is  said  to  be  stationary  in  the  wide  sense  [3.5] . 


We  now  establish  the  following  definition  and  theorem  for  the 


random  processes  z(u)  and  b(u). 


Definition  2: 


The  vector  random  process,  z(u,t),  is  a  two  element  vector  of  the 
random  process,  z(u),  taken  at  time  (u+t)  and  u;  i.e. 


z(  u+t  ) 

z(u, t )  *  ( 3.3 ) 

z(u  ) 


3.1.1.  Invariance  in  the  Mean  with  Time 


Let  z(u)  be  a  zero  mean  (stationary  or  non-stationary )  Gaussian 

random  process,  written  z(u)  =»  z  +  z  (u),  where  z  ,  is  a  random  variable 

o  1  o 

which  we  shall  call  a  reference  phase  and  z  (u)  contains  the  white  and 

1 

independent  increment  (random  walk)  models,  independent  from  z  . 

o 

In  the  treatment  of  each  of  the  phase  random  processes  throughout 
all  other  chapters  of  this  dissertation,  second  order  moments  of  the  rf 
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signal,  b(u),  are  used.  This  requires  two  time  variables  denoted  by  u 

and  u+t.  Ue  also  assume  arbitrary  phase  conditions  at  reference  time,  u, 

whose  variances  are  associated  with  the  multipliers,  c  ,  c  ,  and  c  , 

o  1  3 

for  the  "near  white”  phase,  random  walk  phase,  and  random  walk  frequency 
random  processes  respectively.  It  is  shown  that  these  terms  cancel  out 
in  the  expressions  for  the  auto-correlation  function  of  b(u). 

Since,  in  this  section,  we  are  considering  only  the  first  order 

moment,  we  need  only  one  time  variable  denoted,  u,  with  reference 

time,  0.  For  convenience,  we  include  all  the  terms  constituting  the 

phase  conditions  at  the  reference  time  into  the  term,  z  ,  and  assign  it 

o 

an  arbitrary  variance  multiplier,  c.  The  term,  z  (u),  is  due  to  the 

1 

phase  generated  by  the  independent  white  and  increment  processes 

occurring  from  time,  0,  to  time,  u.  Our  conditions  are; 

z  (0)  -  y(0)  -  z  ,  the  value  of  the  "near  white"  phase  process,  y(u), 

1  o 

at  time,  0,  less  the  arbitrary  phase  component,  z  ,  whose  variance  is 

o 

2  2  2  2 
var(  z  )  »  c  cT  ^c<r+(c-l)d'+c(T. 
o  zoll  o  o32 

Also 

3 

2  2  2  2  2  I  u|  2 

var(z  (u))-<r|u|+or+ccr|ui  +cT  - <r  (u). 

1  1  o  2  2  23  zl 


We  see  that 

2 

var(z  (u))  *  <F  (u)  <  o d,  for  any  |u|  <  oo. 
1  zl 
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Theorea  1:  If  var(z  )  — >  o*>,  then  the  mean  of  the  rf  signal,  b(u) 

o 

is  Invariant  with  time,  l.e. 


j(z(u)) 

E(  b( u  ) )  *  E{  e  } 


TV. 


( 


Proof:  The  mean,  E(  b(u))  is 


E(b(u)) 


1 

\J2  it  <r( u) 


5 


-oc 


exp(  jz 


2 

z 

-  )  dz 

2 

2<r(u) 


E(exp{jz}  ) 


2  2  2 

where  <T(u)^c<J'  +  <T  (u  ). 
z  zo  zl 


In  the  characteristic  function,  F  (p),  if  we  set  p  *  1,  we  obtain 

z 


1 

F  (  1 )  A  E(  exp  {  jz(  u  )}  )  -  exp  { - m  } 

2  “  2  11 


12  12  2 

E(  b(  u  ) )  -  exp{-  -  <r(u)}  *  exp{-  -  [ccT+tr(u)]}  ' 

2  z  2  zo  zl 


Taking  the  following  limit  we  have  the  relationship 


1  i  M 

2  E(  b(  u  )  ) 
c  CT  — >  o*> 

zo 


0 


which  shows  that  E(b(u))  is  invariant  with  u. 


Q.E.D. 


We  know  that  the  oscillator  phase  is  a  renewal  process  and  takes  on 
values  in  the  range  [-IT,  TT  ] .  In  our  discussion  thus  far,  we  have 
assumed  Gaussian  phase  processes  where  the  variable  takes  on  values  in 
the  range  (-ao,  c®  ).  We  should,  therefor s,  restrict  the  phase  variable 
by  the  modulo  2  TT  operation  to  the  range  [-TT,  IT  ] .  We  now  state  the 
following  lemma  which  will  be  obtained  from  Theorem  1. 

Lemma  1 ;  The  average  value 

E(b(u))  =  E(  exp[  jz(u  )]  ) 

is  invariant  with  time  if  the  probability  density  function  of  z  is 

o 

uniformly  distributed  with 

1 

p(z  )  - -  :  |z  |  <  TT  (3.6) 

o  2  T T  o 

Proof:  From  theorem  1,  one  needs  only  to  convert  the  variable  of  the 
Gaussian' density  function  to  modulo  2  TT. 

For  a  particular  point,  say  x  within  -TT  <  x  <  TT,  the  modulo 

2  TT  value  is  the  superposition  of  all  z  =■  x  +  n  2IT;  where  n  is 

o 
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n—cffi,  ...-2,  -1,  0,  1,  ...on,  as  shown  in  figure  3.1.  In  terms  of  the 
original  Gaussian  function, 


P(z  ) 

0 


1 


n  c  <r 
zo 


(x  -  n2IT) 


2 

2  c  cr 
zo 


I 

I 

J 


Multiplying  and  dividing  through  by  <T 

20 


I  ITT  we  have 


p(z  ) 

o 


1 

2  IT 


CK> 


ZO  1 

_____ —  y  exp{ — 

\T3tt  <r  2c 


zo  n—oc 


2ITn 


zo 


l 

Then  from  the  approximating  function  and  setting 


2  IT 


2  IT 


Au  and  - n  «  u,  we  have  for  I  z  I  <  IT _ 

[T  <r  o 


\J  zo 


j  |  M  p(  z  ) 

2  o 

C<T  . >OP 

zo 

/\u — -K) 


1^ 

The  approximating  function  approach  allows  integration  by  the 
summation  of  discrete  values  in  the  limit.  The  details  of  the  approach 
can  be  found  in  [3.6]. 
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FIGURE  3.1  SUPERPOSITION  OF  pdf  SEGMENTS  OF  2  IT  FROM 
GAUSSIAN  DENSITY  FUNCTION 

* 


It  should  be  noted  that  an  oscillator  which  has  been  drifting  over 
a  long  period  of  time  will  have  a  reference  phase  which  is  completely 
uncertain,  i.e.  it  has  a  uniform  distribution  over  2  IT  radians.  This 

i 
) 

I  should  be  an  obvious  condition  that  one  would  expect. 
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We  can  also  find  the  average  value  of  b(u)  by  using  the  density 
function  of  equation  (3.6)  as  follows 


rr 


E(b(u))  -  ^ 


exp{jz(u  )} - dz(u  ) 

2  rr 


1 

j2  rr 


-rr 


rr 

exp  { jz } 

-rr 


exp{jIT}  -  exp{-jTT}  2  sin(IT) 


(3.8) 


J2  rr 


2  rr 


which  agrees  with  the  result  fron  equation  (3.5). 


3.1.2.  Invariance  in  Auto-correlation  Function  with  Time,  u 

Let  the  vector  random  process  z(u,t)  ■  of  real  elements  have 
a  covariance  matrix  of  the  form 


T 

E(  z  z  ) 


a(  u  )  +  g(  t  ) 
a(u)  -  h(t  ) 


a(u)  -  h(t  ) 
a(u ) 


(3.S) 


The  rf  oscillator  signal,  b(u),  is  given  by 

b(u)  »  exp{  j  z(u)  }  (3.10) 

where  z(u)  is  a  Gaussian  phase  random  process  which  is  not  necessarily 
stationary  and  u  denotes  real  time.  The  auto-correlation  function  is 
given  by 


R  (t)  -  E(  exp {  j[  z(u+t)  -  z(u)  ]}  ) 


(3.11) 


Theorem  2;  The  auto-correlation  function  of  the  rf  oscillator  signal, 
b(u),  with  the  covariance  matrix  given  by  equation  (3.9)  depends 
only  on  t  for  all  u  and  is  given  by 

*  g(t) 

R  (t)  =  E(  b(u+t)  b  (u)  )  *  exp{  - - -  h(t)  }  (3.12) 

b  2 


Proof:  We  proceed  to  find  the  auto-correlation  function  of  the  waveform, 
b(u),  by  using  the  characteristic  function  approach  of  Appendix  A. 

When  b(u)  =  exp  {  jz(u)},  with  z_  Gaussian, 

* 

R  (t)  A  E{  b(  u+t  )  b  (u)  } 
b 


R  (t)  -  F  (  +1,  -1  )  (3.13) 

b  z 


where  F  is  the  characteristic  function  of  z_,  i.e.  for  a  Gaussian  random 
vector,  z_,  with  covariance  matrix,  M 

z_ 

1  T 

F  (jO  -  exp { - £  M  £  } 

z  2  z_ 

T 

If  The  covariance  matrix  elements  are  m  ,  then,  for  p  ■  [  1,  -1  ] 

ij 

1 

R  (  t  )  ■  exp  {  -  -  [  m  -  2m  +  m  J  } 

b  2  11  12  22 


When 


M 

z 


a(u)  +  g(t) 
a(u  )  -  h(t) 


a(u  )  -  h(  t  ) 
a(u  ) 


we  have  finally 


(3.14) 


-  55  - 


g(t ) 

R  (  t  )  *  exp  {- - -  h(  t  )} 

b  2 


(3.15) 


Note  that  although  z(u)  is  non-stationary ,  the  auto— correlation  function 

of  the  oscillator  rf  signal,  b(u),  depends  only  on  t  for  all  u. 

Q.E.D. 

Next,  we  investigate  the  ergodic  properties  of  our  phase  models. 

3.2.  Ergodicity: 

It  is  well  known  [3.1],  [3.2]  that  if  a  process  is  ergodic,  then 
time  averages  are  equal  to  the  corresponding  ensemble  averages.  The 
properties  of  the  process  can  therefore  be  determined  from  one  sample 
function  which  is  useful  from  the  practical  standpoint. 

In  this  section  we  establish  a  number  of  theorems  on  ergodicity  for 
three  statistical  parameters  of  the  rf  oscillator  signal ;  the  mean,  the 
auto— correlation  function,  and  the  power  spectral  density.  We  consider 
all  the  three  phase  random  processes; the  white  phase,  random  walk  phase, 
and  the  random  walk  frequency. 


3.2.1.  Ergodicity  in  the  Mean 


A  test  for  ergodicity  in  the  mean  of  a  random  process  is  stated  as 
follows  [3.2],  pp  328 
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Theorem: 


lim 

T  ■— ’frcg) 


1 

2  T 


T 


b(  u  )  du 


E(  b(  u  )  )  = 


(3.16) 


Iff 

2T 

1  /  lt|  2 

lim  -  \  (  1 - )[  R  (t)  -  ^  ]  dt  =  0  (3.17) 

T> — »qp  T  /  2T  b 

0 


* 

where  R  (t)  =  E[b(u+t)b  (u)]  is  the  auto-correlation  function  of  b(u). 
b 

The  mean,  and  the  auto-correlation  function,  Rb( t ),  of  the  phase 
random  processes  have  been  determined  in  Appendix  D.  These  are  required 
for  the  invest!  ation  of  ergodicity  in  the  mean.  First  we  consider  the 
white  (or  "near  white”)  phase  process. 

Theorem  3: 

The  rf  signal  process  due  to  "near  white"  phase  instability  is 
ergodic  in  the  mean  if  the  random  phase  at  any  time,  u,  approaches 
complete  decorrelation  after  some  time  shift,  t,  i.e.  for 

0  <  T  /  2  <  |  1 1  <c*> 
o 

the  auto-correlation  function  satisfies 

* 

R  (t  )  A  E(y(u+t>/  (u))  =  o  (3.18) 

y 


Proof:  Using  appendix  A  and  appendix  B  it  can  be  shown  from  the 

characteristic  function  that  for  the  white  phase  process,  y(u),  with 


I  r  I  > 


t  /  2, 

o 


the  auto-correlation  of  the  rf  signal,  b(u),  is 


1 

E(  exp  {  i[y(u+t  )-y(u  )J  }  )  =  F  (  1,-1  )  =  exp{  -  -(m  +  m  )+m  } 

y  2  11  22  12 


2 

=  R  (  t  )  =  exp(  -  o  ) 
b  o 


T 

o 

t  > -  (3.19a) 

2 


Also  from  the  characteristic  function  approach  it  is  easy  to  show  that 


1 

F  ( 1 )  =  exp {  -  -  m  } 
y  2  11 


2  o 

7\_  =  E(  exp(  jy(u  ))  *  exp(-  <r - )  (3.19b) 

b  o  2 


Substitution  of  equation  (3.19)  into  (3.17)  gives 


2T 

1  (  t  2  2 

lim  -  \(1  - - )[exp(-  <?  )  -  exp(-  <T  c  )]  dt  =  0 

T — T  /  2T  o  o  o 

0 


c  =  1 
o 

(3.20) 


2  2 

=  exp(  -  <T  )  -  exp(  -  cr  c  )■»*»  0  :  c  1 

o  o  o  o 


i.e.,  the  rf  signal  process  is  ergodic  in  the  mean  if  for  some  t  >  T  /2, 

o 

c  -  1  ■  0.  The  covariance  matrix  for  this  condition  is 
o 


T  2 

E(  yy  )  =  cr 


T  /2  <  I  tl  <  op  (3.21) 
o 


Since  the  off  diagonal  terms  denote  correlation  with  respect  to  the  time 
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shift,  t,  then  c  =  1  implies  that  the  phase  process,  y(u),  has  become 

o 

decorrelated  at  time  u+t.  Hence, 

c-1— *R(t)«0  :  |  t|  >  T  /2. 

o  y  o  Q.E.D. 

If  a  correlated  component  of  the  phase  at  time,  u,  does  exist,  it  must 
be  known  or  deterministic  in  order  for  ergodicity  of  the  mean  to  apply. 

Theorem  4: 

The  rf  signal  process,  b(u),  due  to  the  random  walk  phase 
instability  process  is  ergodic  in  the  mean  if  the  reference  phase  at 
time,  u,  is  uniformly  distributed  phase  over  2IT;  c. f. ,  Theorem  1, 
Lemma  1 

Proof:  It  has  been  shown  in  the  appendices  that  for  the  random  walk 

phase  process,  x(u), 

2  t 

R  (  t  )  =  exp(-  <r - )  t  >  0  (3.22a) 

b  1  2 

Using  the  characteristic  function  in  conjunction  with  the  arbitrary 
phase,  we  can  evaluate 

2 

ft  =  exp(-  (Tc/2)  (3.22b) 

b  1  1 

Substituting  equation  (3.22)  into  equation  (3.17)  gives 


1  (  t  2  t  2 

lim  -  \  (1 - )[exp(-  cr  - )  -  exp(-  cr  c  )]  dt 

T*— >a*>  T  )  2T  12  11 

0 


2 

-  exp(-  <r  c  ) 
1  1 


0  :  c  <  a» 

1 

■  0  :  c  =  cm 

1 


(3.23) 


This  completes  the  proof  for  ths  random  walk  phase.  Q.E.D. 


Theorem  5: 

The  rf  signal,  b(u),  due  to  the  random  walk  frequency  process,  v(u) , 

Is  ergodlc  In  the  mean  if  the  reference  phase  at  time,  u,  is  uniformly 
distributed  over  2  IT,  cf.  lemma  1. 

• 

Proof :  For  the  "Random  Walk"  Frequency,  v(u),  it  can  be  shown  from  the 
appendices  that  for  the  process,  v(u), 

2  2 

R  (t)  -  exp(-  <r  [3c  +  t]  t  /6)  (3.24* \ 

b  2  2 

Using  the  characteristic  function  in  conjunction  with  the  arbitrary 
phase  variance  we  find  the  mean  of  the  rf  signal 

2 

X  *  exp(-  <r  c  / 2)  (3.24b) 

b  2  3 


Substituting  equation  (3.24)  into  equation  (3.17)  gives 


rrn 


lim 
T— — >oo 


t  2  2  2 

—  )[exp(-  c r  [3c  +  t]t  /6  )  -  exp(-  <r  c  /2)]  dt 
2T  2  2  2  3 


2 

=  -  exp(-  cr  c  /2)  0  :  c  <  c® 

2  3  3 


=0  :  c  =  oo  (3.25) 

3 

Q.E.D. 

This  completes  the  proof  for  the  random  walk  frequency  process,  v(u). 

We  now  turn  our  attention  to  ergodicity  of  the  auto-correlation 
function. 


3.2.2.  Ergodicity  in  Auto-correlation  Function 


The  three  phase  instability  models,  white  phase,  random  walk  phase, 
and  random  walk  frequency  will  now  be  considered  for  the  ergodicity  in 
auto-correlation  function  of  the  oscillator  rf  signal. 

A  test  for  ergodicity  in  the  auto-correlation  function  of  a  random 
process  is  stated  as  follows  [3.2]  (pp  330). 


Theorem:  For  a  given  s 


lim 
T — >o© 


iff 


- 

2T  ) 


b(u+s)b  (u)  du  =  E{b(u+s)b  (u)}  =  R  (s) 

b 


(3.26) 


-T 


2T 


lim 

T— — >a© 


1  f  ' 

-  \  (  1 - )[R 

T  /  2T  i 


(t)  -  R  (s)]  dt  =  0 
#  b 


(3.27) 
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* 

where  R  (s)  =  E[b(u+s)b  (u)j  is  the  auto-correlation  function  and 
fc 

*  * 

R  (  t  )  =  E{  b(u+s+t)b  (u+t)b  (u+s)b(u)  ]  is  the  fourth  order  moment  of 

H 

the  rf  signal,  b(u).  We  will  proceed  to  establish  theorems  relating  to 
the  ergodicity  of  the  auto-correlation  function  of  the  process 

b(u  )  *=  exp(  jz(u  )). 

To  accomplish  this  we  need  to  determine  the  second  and  the  fourth  order 
moments  appearing  in  equation  (3.27).  We  will  use  the  procedure  and 
results  given  in  appendix  D  to  determine  these  moments. 

Theorem  6; 

If  the  phase  random  process  of  an  rf  signal,  b(u),'  is  a  white  noise 
phase  process,  y(u),  then  b(u )  is  ergodlc  in  its  auto-correlation 
function. 

Proof:  With  appropriate  substitutions  of  R  (t)  and  R  (s)  from 

b 

Appendix  D,  for  the  white  phase  process,  equation  (3.27)  is  evaluated 
under  three  regions  for  the  shift  variables,  t  and  s  i.e. 

1)  0<t  <8+T/2  sails 

o 

and  2)  0<s  +  T/  2<t  :s<T/2 

o  o 

3)  0<s  +  T/  2<t  :  s  >  T  /  2 

o  o 

The  integral,  equation  (3.27),  for  the  first  region  is  finite  due  to  the 
fact  that  the  area  under  the  white  (or  "near  white")  impulse  type 
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function  is  finite  for  a  finite  interval,  T  . 

o 


For  the  second  region,  we  have 


2T 


lim 

T - >003 


-  (  (1  --- 
T  )  2T 


2  4  2  4 

)[  exp(  -  O' - s  )  -  exp(  -  <T - s  )]  dt 

o  T  o  T 

o  o 


2T 


lim 
T — fcm 


-  ^  Cl  -  — 

T  /  2T 


)[0]  dt  -  0 


(3.28a) 


:  s  <  T  /2  <  t-s 
o 


and  for  the  third  region,  we  have 


2T 


lim 
T — >ao  T 


2  2 

- )[  exp(  -  2 <r  )  -  exp(  -  2o~  )  ]  dt 

2T  o  o 


2T 


lim 
T—  )oc 


T  )  2 


)f0]  dt  -  0 


(3.28b) 


2T 


:  T  /2  <  s  <  t-T  /2 
o  o 


We  see  that  the  limit  of  equation  (3.27)  approaches  zero  for  the  three 
conditions  above.  Hence,  equation  (3.27)  is  satisfied  for  all  s. 

Q.G  *D . 

We  conclude  that  the  auto-correlation  function  of  the  rf  signal,  b(u), 
is  ergodic  when  the  phase  instability  is  a  white  (or  near  white)  random 
process,  y(u). 


Theorem  7: 

If  the  phase  random  process  is  a  random  walk  pru..ess,  x( u ),  then 
the  rf  signal,  b(u),  is  ergodic  in  auto-correlation  function. 


Proof:  The  appropriate  second  and  fourth  order  moments  for  the  random 
walk  phase  are  taken  from  appendix  D  and  are  substituted  into 
equation  ( 3.27 ). 

It  can  be  argued  when  t  is  such  that  0  <  t  <  s,  the  integral  of 
equation  (3.27)  is  finite  and  so  equation  (3.27)  approaches  zero.  For 
all  other  values  of  t  such  that  0  <  s  <  t,  the  integral  is  identically 
zero.  That  is 

2T 

lft  2  2 

lim  — « •  \  (  1 - )[exp(-  o' s  )  -  exp(-  0's)]  dt 

T— T  /  2T  1  1 

0 

2T 

1  (  * 

*  lim  —  V  (1 - )[0]  dt  -  0  (3.29) 

T — T  J  2T 

Therefore,  equation  (3.27)  is  satisfied.  Q.E.D. 

We  conclude  that  the  the  rf  signal,  b(u),  is  ergodic  in  auto-correlation 
function  of  the  when  the  phase  random  process  is  a  random  walk  phase 
process,  x(u). 

Theorem  8: 

If  the  phase  random  process  is  a  random  walk  frequency  process , 
v(u ),  then  the  rf  signal,  b(u ),  is  not  ergodic  in  auto-correlation 
function. 

Proof:  For  random  walk  (Brownian  Motion)  frequency  instabil  y, 
equation  (3.27)  uses  the  fourth  and  second  order  moments  obtained  in 
equation  (D.7).  The  result  is 


2T 


lim 

T — >op  T 


2  2 
t  2  s  2  s 

■  — *  )[exp(-  tr  —  [  3 1— s ]  )  -  exp(-  <T  —  [3c  +s]  )J  dt 
2T  2  6  2  3  2 


The  first  term  in  the  brackets  has  a  negative  exponent  and  will  bring 
the  limit  to  zero.  The  second  term  is  a  constant  which  leads  to 


1 


2T 


lim 

T — >o© 


T  J  2T 


2  s 

)[-exp(-  <T  —  [3c  +s]  )]  dt 
2  3  2 


2  s 

-exp(-  O’  —  [3c  +s]  0 

2  3  2 


(3.30) 


We  conclude  that  the  rf  signal,  b(u),  is  not  ergodic  in  auto-correlation 

function  of  the  when  the  phase  random  process  is  a 

random  walk  frequency  process,  v(u).  Q.E.D. 

Observe  the  rf  signal,  b(u),  is  not  ergodic  in  auto-correlation 

function  when  the  phase  random  process  consists  only  of  a  long  term 

drift  offset  frequency.  This  is  obvious  due  to  the  presence  of  c  and 
2  2 
<T*  in  equation  (3.30). 

2 

For  certain  processes,  note  that  ergodicity  may  exist  in  the 
auto-correlation  function  and  yet  may  not  exist  in  the  mean.  This  is 
easily  seen  through  the  fact  that  the  mean  phase,  though  random  and 
non-ergodic  across  the  ensemble,  is  subtracted  out  for  any  one  sample 
function  in  the  process  of  finding  the  auto-correlation  function.  This 
is  done  through  the  functional  operation 


R  (  t  )  *  E(  b(u+t  )b  (u  ))  -  E(  exp[  jz(u+t  )-jz(u  )j  ) 
b 
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eliminating  any  ef feet, a  which  the  mean  phase  variable  has  on  the  rf 
auto-correlation  function. 

This  effect  is  brought  out  in  the  next  theorem  on  ergodicity  of  the 
power  spectral  density.  In  this  case  the  specular  component  and  spread 
component  are  not  necessarily  mutually  ergodic. 

3.2.3.  Ergodicity  in  Power  Spectral  Density 

Many  references  [3.2],  [3.7]  provide  conditions  for  determining  the 
power  spectrum  from  the  average  random  power, i.e.  the  time  average. 
When  these  conditions  are  met,  then  the  limiting  value  of  the  time 
average  is  equal  to  the  power  spectral  density,  i.e. 

1 

S(  w )  =*  lim  - 

T  — >  op  2T 

and  the  power  spectral  density  is  ergodic. 

Ergodicity  of  the  power  spectral  density, 

op 

S(w)A  ^  R(  t  )  exp(-jwt)  dt 
— GtD 

requires  first,  that  the  expected  value  of  the  "average  random  power”, 

T  2 

1  I  ( 

S  (w)  ^  -  \  b(u)  exp(-jwu)  du 

T  2T  I  / 

-T 


T  2 

x(u)  exp(-jwu)  du 
-T 


should  tend  to  S(w)  in  the  limit;  and  second,  that  its  variance  should 


tend  toward  zero  in  the  limit,  T  — >  a®.  We  will  only  address  the  first 
of  these  requirements  which  will  provide  us  with  a  necessary  condition 
for  testing  ergodicity  in  power  spectral  density.  We  proceed,  using  the 
theorem  provided  on  page  343  of  [3.1]  which  follows 


Theorem:  If 


5 


-a© 


I  tR  (  t  )|  dt  <  c© 
b 


then 


lira  E{S  (w)} 
T — >o®  T 


“  S(w) 


a® 


■5 

-0© 


R  ( t )  exp(  -jwt )  dt 
b 


T  2 

where  1  I  f 

S  (w)  - -  \  x(u)  exp(-jwu)  du 

T  2T  I  / 

-T 


(3.31) 


(3.32) 


* 

and  R(  t  )  »  E{x(u+t)x  (u)}. 

That  is,  if  equation  (3.31)  is  satisfied,  then  a  necessary  condition  has 
been  net  for  the  power  spectral  density  of  the  process,  b(u),  to  be 
equal  to  the  time  average  of  the  power  spectrum  obtained  from  absolute 
value  squared  of  the  Fourier  transform  of  the  sample  function.  We  now 
state  our  theorem  for  ergodicity  of  Power  Spectral  Density. 
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Theorem  9: 


If  a  random  phase  process,  z(u),  consists  of  the  sum  of  a  white 
Gaussian,  finite  variance ,  phase  process  and  higher  order  independent 
increment  processes ,  then  the  corresponding  power  spectral  density  is 
not  ergodic.  The  random  correlated  component  exists  over  the  interval, 
-<xD  <  u  <  oo ,  and  has  a  finite  Gaussian  variance,  fl"(white),  i.e.  not 
uniformly  distributed  over  modulo  2TT,  c. f.  Lemma  1. 


Proof:  It  has  been  established,  in  our  work  thus  far  and  elseware 
[3.3],  [3.8],  [3.9]  that  the  auto-correlation  function  of  an  oscillator 
rf  signal  containing  white,  random  walk  phase,  and  random  walk  frequency 
components  is  of  the  form 


2  2  2  3 

R  (  t  )  =  exp(-<f  —  Itl  —  b|t|  —  c|tl  —  d|t|  )  :  |  t|  _<  T  /2 

b  o  To  o 

(3.33) 

2  2  3 

R  (  t  )  =  exp(-<T  -  b|t|-c|t|-d|ti  )  :  |t|>T/2 

bo  o 


Clearly,  the  exponential  terms  containing  |t|  and  powers  there-of 
force  the  integral  in  equation  (3.31)  to  converge  to  some  finite  value 
so  that  equation  (3.32)  is  satisfied.  The  only  condition  we  will 
investigate  is  when  b  =  c=  d*'0,  i.e.  the  random  walk  phase,  the 
offset  frequency,  and  the  random  walk  frequency  component  are  absent  (or 
of  zero  level).  Then  only  the  white  phase  component  is  present.  We 
have  for  this  case 
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Defining 


2  2 

R  (  t  )  -  exp(  -«T  —  I  1 1  ) 

b  o  T 

o 


2 

R  (  t  )  -  exp(  -<r  ) 
b  o 


I  t|  <  T  12 
o 


I  1 1  >  T  /2 
o 


2  2 

f(t  )  *  rect{t/T  }  [  exp( -<T - |t|) 

o  o  T 

o 

where  rect{t}  M  :  I  t|  <  1/2 

^  0  :  otherwise 

we  can  write  the  above  relationships  as  follows: 

2 

R  (t)  -  f(t)  +  exp(  -<r  )  (3.34) 

b  O 

Substituting  equation  (3.34)  into  equation  (3.31),  we  have 

dt 

(3.35) 

Since  f( t )  is  bounded  and  is  non-zero  for  finite  duration  only,  the 
first  integral  on  the  right  of  equation  (3.35)  is  finite  and  therefore 
satisfies  the  necessary  condition  stated  in  equation  (3.31).  It  remains 
to  determine  the  effects  which  the  second  integral  on  the  right  of 
equation  (3.35)  have  on  the  ergodicity  of  the  power  spectral  density. 

This  integral  arises  fror  that  tern  of  the  "near  white"  phase 
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o© 
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I  2 

1  t[f(t  )  +  exp(-<T  )] 

dt  ■  \ 

t  f(t) 

C  1  2 

dt  +  \  It  exp(-<T  ) 

1  o 
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-  exp(  -<T  )  ] 
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[T. 


re 


re 


2 

component,  equation  (3.34),  which  contributes  the  constant,  exp{-C*}» 

o 

to  the  auto-correlation  function  for  all  T  /2  <  | tl  <  ao.  Following 

o 

through  with  the  evaluation  of  the  second  integral,  in  equation  (3.35) 

2 

(equation  (10-27)  of  [3.2]  )  for  the  pedestal,  exp(-<T  ),  we  obtain 

o 

2  1  -  cos(2Tw) 

lim  E[  S  ( w )]  “  lim  exp(  -  0"  ) - 

T— >qo  T  T  ■  — >op  o  2 

T  w 

As  T— frxp,  this  becomes  an  impulse  function  in  w,  ([3.9],  pp  378, 
equation  633).  We  can  write 

2 

lim  [S  (w)]  -  2TT  exp(  -  o'  )  d(w)  (3.36) 

T«-— >ap  T  o 

where  d(w)  is  an  Impulse  function. 

We  see  how  the  pedestal  component,  i.e.  the  constant 
auto-correlation  function  for  all  t,  of  the  "near  white”  phase  process 
gives  rise  to  the  specular  component,  or  the  frequency  component 
representing  the  correlated  random  part  of  the  oscillator.  The  pedestal 
component  is  not  ergodic  in  power  spectral  density.  The  spectral 
spreading  due  to  the  time  decorrelated  component,  the  first  term  in 
equation  (3.35),  does  however  satisfy  the  condition  for  ergodicity  given 
in  the  above  theorem.  Hence,  theorem  9  is  proved.  These  properties  are 
shown  in  figure  3.2. 

Theorems  have  been  provided  for  stationarity  and  ergodicity  of  the 
rf  signal,  b(u),  when  the  phase  random  process  is  white,  random  walk, 
and/or  when  the  frequency  random  process  is  random  walk.  Properties 
covered  were  the  mean,  auto-correlation  function,  and  power  spectral 
density.  Table  3.1  summarizes  the  results. 
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FIGURE  3.2  WHITE  PHASE  INSTABILITY  PROPERTIES 


TABLE  3.1  CONDITIONS  FOR  ERGODICITY 
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MEAN 
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ERGODIC: 

RANDOM  COMPONENT 
OF  STARTING  PHASE 
TENDS  TOWARD  COMP¬ 
LETE  DECORRELATION 
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HAVING  ANY  pdf 


AT  u  <  oo 
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x(u) 

RANDOM 
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PHASE 
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STARTING  PHASE  HAS 
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v(u) 

RANDOM 
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POSSIBLY  ERGODIC: 
ALL  CONDITIONS 


POSSIBLY  ERGODIC 
ALL  CONDITIONS 


REFERENCES 


[3.1]  Davenport,  W.  B.  and  Root,  W.  L.,  "Random  Signals  and  Noise", 
McGraw-Hill.  New  York,  1958,  sec  6-6  pp  107-108,  sec  12-3 

pp  257-259. 

[3.2]  Papoulis,  A.,  "Probabiility ,  Random  Variables  and  Stochastic 
Processes",  McGraw-Hill  1965,  pp  328,  330,  343-344. 

[3.3]  Rutman,  J.,  "Charactorization  of  Phase  and  Frequency  Instabili¬ 
ties  in  Prescision  Frequency  Sources:  Fifteen  Years  of  Progress", 
Proc.  IEEE,  vol  66,  no  9,  pp  1048-1075,  September  1978. 

[3.4]  Boileau,  E.  and  Picinbono,  B.,  "Statistical  Study  of  Phase 
Fluctuations  and  Oscillator  Stability”,  IEEE  Trans,  on 
Instrumentation  and  Measurements,  Vol  25,  No.  1,  Mar  1976, 
pp  66-75. 

[3.5]  Davenport,  W.B.,  Jr.,  "Probability  and  Random  Processes", 
McGraw-Hill,  New  York,  1970,  Chapter  9,  Section  7. 

[3.6]  Davenport,  W.B.,  Jr.,  "Probability  and  Random  Processes", 
McGraw-Hill,  New  York,  1970,  Chapter  7,  Section  5. 

[3.7]  Sage,  A.P.,  Melsa,  J.L.,  "Estimation  Theory  with  Applications  to 
Communications  and  Control",  McGraw-Hill,  New  York,  Sec.  3.4,  1971. 

[3.8]  Barnes,  J.  A.,  et  al,  "Characterization  of  Frequency  Stability", 
IEEE  Trans.  Instrum.  Meas.,  vol.  IM-20,pp.  105-120,  May  1971. 

[3.9]  CRC  Standard  Math  Tables,  25th  Edition,  CRC  Press,  West  Palm 
Beach,  Florida,  1978. 


-  73  - 


IV  SPECTRAL  SPREADING  DUE  TO  PHASE  INSTABILITY 


4.0 .  Introduction 

It  is  well  known  that  a  monochromatic  oscillator  signal  provides  a 
specular  ( impulse )  mapping  into  the  spectrum  domain.  Such  a  signal, 
when  modulated  by  an  envelope  waveform,  fo(u),  provides  a  frequency 
domain  convolution  of  the  specular  spectrum  with  the  envelope  sprectrum 
to  result  in  the  modulated  system  spectrum.  Normally,  the  ideal 
envelope  spectrum  is  characterized  by  true  nulls,  of  infinite  depth,  in 
the  frequency  domain.  When  convolved  with  a  frequency  specular  the 
spectrum  of  the  modulated  system  waveform  reproduces  these  infinitely 
deep  nulls.  With  such  an  ideal  spectrum  the  potential  exists, 
theoretically;  to  process  a  second  signal  located  at  one  of  these  nulls 
to  allow  almost  perfect  discrimination  between  the  two  signals  when  only 
low  level  thermal  noise  is  present.  A  simple  example  of  this  situation 
is  found  in  radar  systems.  When  the  transmitted  signal  of  a  ground 
based  radar  is  such  a  modulated  waveform  the  return  from  a  moving  target 
is  doppler  shifted  relative  to  the  return  from  terrain  clutter.  By 
observing  the  return  within  the  nulls  of  the  unshifted  spectrum  one  can 
detect  the  moving  target  against  clutter  with  a  level  of  performance 
limited  only  by  the  thermal  noise  of  the  system. 

With  the  introduction  of  phase  instability  the  spectrum  of  the 
oscillator  is  no  longer  a  specular  impulse.  The  specular  broadening 
effect  due  to  the  oscillator  instability  causes  a  spectral  spreading  of 
the  system  modulated  waveform  with  ar  associated  deterioration  of  the 
null  depths.  As  the  potential  discrimination  in  the  null  region 
degrades  due  to  phase  instability  whi-h  is  in  addition  to  thermal  noise, 
the  system  performance  suffers  further.  As  mentioned  earlier,  we  only 
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consider  the  effect  of  phase  instability  in  this  dissertation.  The 
degradation  is  interrelated  with  a  variety  of  system  parameters  such  as 
the  observed  spectral  region;  the  level  and  type  of  instability,  i.e. 
white,  Wiener,  etc.;  and  the  type  of  envelope  waveform  used  to  modulate 
the  oscillator.  We  will  discuss  the  interrelation  between  these 
parameters  and  their  spectral  spreading  effects  which  deteriorate  the 
observable  nulls. 

Nest,  we  will  describe  the  modulation  system  and  the  operations 
involved.  The  modulated  waveform,  denoted  by  fk(u)  in  equation  (1.2), 
(henceforth  denoted  by  f(u))  with  average  autocorrelation  function 
Rf(  t  )  in  equation  (1.4),  (henceforth  denoted  by  R(t))  is  often  perceived 
s  a  train  of  pulses  as  shown  in  figure  4.1.  The  waveform,  f(u),  is  the 
output  of  the  modulator  shown  in  figure  1.1a  and  its  average 
auto-correlation  function,  R(  t  )  is  produced  by  the  equivalent  diagram  of 
figure  1.1b.  The  waveform,  f(u),  consists  of  a  stationary  random 
process,  b(u  ),  driven  by  the  non-stationary  random  phase  process,  z(u), 
shown  in  figure  4.2a. 


— *4  k"  P 

JL. 

N  -I 
2 


FIGURE  4.1  ASSUMED  PULSE  TRAIN 
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The -process  is  multiplied  by  some  ideal  deterministic  waveform,  for 
example,  a  finite  train  of  N  narrow  pulses,  fo(u),  shown  in  figure  4.2b. 
This  is  generally  a  non-stationary  random  process  which  produces  the 
pulse  train,  f(u),  whose  amplitude  and  random  phase  is  shown  in 
figure  4.2c.  Finally  this  pulse  train  is  convolved  with  the  pulse 
shape,  so(u),  shown  in  figure  4. 2d  to  produce  the  desired  train,  s(u), 
shown  in  figure  4.2e.  If  the  pulse  shape  were  rectangular  it  would 
produce  the  pulse  train  shown  in  figure  4.1.  Once  again  in  this 
chapter,  we  shall  treat  white  phase  and  random  walk  phase  models 
seperately  from  the  random  walk  frequency  model.  Only  analytical 
results  are  obtained.  Numerical  results  will  be  treated  in  the  next 
chapter. 

4.1.  White  Phase  and  Random  Walk  Phase 

As  indicated  earlier,  for  this  model  the  oscillator  signal,  b(u), 
includes  a  non-stationary  phase  random  process,  z(u),  which  is  the  sum 
of  a  "near  white"  process,  y(u),  and  an  independent  Wiener  (BM  or  random 
walk)  process,  x(u). 

In  this  section,  we  shall  treat  the  following  three  modulating 
envelope  waveforms: 

1)  A  cw  waveform  as  given  in  equation  (2.5a); 

2)  An  infinite  train  of  narrow  pulses,  d(u-iT); 

3)  A  finite  train  of  N  short  pulses  as  given  in  equation  (2.5b). 


These  waveforms  are  shown  in  figures  4.3a,  4.3b,  and  4.3c 

respectively . 
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FIGURE 


FIGURE  4.2.  MODULATION  OF  AN  OSCILLATOR  RF  SIGNAL 
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<« 


A 

1  c 

R(t)-lim  -  \  f(u)f(u-t)du  (4.1a) 

o  A— >od  2A  )  o  o 

-A 


For  the  finite  pulse  train  of  short  pulses,  we  use  the  time  average 
auto-correlation  function  defined  as: 

o£> 

7  * 

R  (  t  )  =  \  f  (u  )  f  (u-t )  du  (4.1b) 

o  /  o 

-go 

Multiplying  equation  (4.1a)  and  equation  (4.1b)  by  the  auto-correlation 

function  of  the  oscillator  rf  signal,  R  (  t  ),  and  taking  the  Fourier 

b 

transform  of  the  product  results  in  the  average  power  spectrum  and 

average  energy  spectrum  respectively  [4.1].  The  average  power  spectrum 

becomes  power  spectral  density  only  when  the  modulating  envelope 

waveform,  f  (u),  is  stationary,  e.g.  a  cw. 
o 

From  equation  (1.4),  we  can  obtain  the  average  auto-correlation 
function  of  our  system  output  wav  form, 

R(t)  -  R  (t)  R  (t)  (4.2). 

b  o 

Recall  that  the  system  output  waveform  is  generated  by  taking  the 
following  product,  equation  (1.2) 

f(u)  =  f  (u  )  b(u  )  (4.3). 

o 

The  general  procedure  that  we  shall  follow  for  each  of  the  three 
waveforms  is: 

i)  Evaluate,  the  oscillator  auto-correlation  function,  Rb( t ), 
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using  the  procedures  developed  in  chapter  two,  equations 
(2.10)  and  (2.15). 

ii )  Evaluate  the  time  average  auto-correlation  function,  Ro(t),  of 
the  modulating  envelope  waveform,  fo(u),  using  equation  (4.1). 

iii )  Determine  the  average  auto-correlation  function,  R(t)  of  the 
system  output  using  equation  (4.2)  by  substituting  Rb( t )  and 
Ro(  t  )  found  in  (  i  )  and  (ii)  above. 

iv)  Find  the  average  power/energy  spectrum  through  the  Fourier 
transform  [■'.!],  [4.2]. 


-jwt 


S(w) 


R(t) 


dt 


4.1.1.  The  cw  envelope  waveform: 

This  waveform,  f  (u)  =  c,  is  a  stationary  random  process.  The 
o 

av  rage  auto— correlation  function  of  the  output  waveform,  f(u),  becomes 
the  auto-correlation  function  and  the  average  power  spectrum  becomes  the 
power  spectral  density  [4.1] .  For  the  cw  envelope  waveform,  the  time 
average  auto-correlation  function  is  the  auto-correlation  function 
(ensenfele  av  rage)  and  is  given  by: 


1 

R  ( t  )  =  lim  - 

o  A— >o©  2A 


c  c 


dt 


(4.4). 


Using  equation  (4.2)  we  determine  the  auto-correlation  function  of  the 
modulated  output  waveform  as  follows: 
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where 
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The  shape  of  this  function  is  shown  in  figure  2.2. 

A  generalization  of  the  above  equation  can  be  made  to  include  other 
shapes  for  the  Impulse  type  auto-correlation  function.  By  simply 
introducing  a  superposition  of  any  nuofcer,  say  (L-l),  of  "near  white" 
phase  processes,  each  having  different  triangular  slopes  and 
decorrelation  times,  see  figure  4.4,  we  can  extend  the  stationary  "near 
white"  phase  model,  y(u),  to  include  a  variety  of  auto-correlation 
functions  consisting  of  concave  curvatures.  It  is  not  difficult  to  show 
that  the  expression  for  these  extended  cases  can  be  written  as 


R(t)  ■  | cl  exp{-(k'  +  k  !tl  )} 

i-1  i 


(4.5b) 
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where  i  -  0,  1 ,  2, 


•  •  •  •  ( L—  1 ) »  L ; 


T  -0<T<T*oo,  and  for 
-1  i  L 


T  T 

i-1  i 

-  <  lt|  <  -  , 

2  “2 


2 


T 


j 


and  Che  T's  are  the  decorrelation  times  of  the  additive  “near  white” 
i 

Gaussian  phase  random  processes.  The  logarithm  of  equation  (4.5b)  is 
plotted  in  figure  4.4. 


In  Rf(t)**CKj.,*Kj  Itl) 


FIGURE  4.4.  LOGARITHM  -  OSCILLATOR  PHASE  AUTO- CORRELATION  FUNCTION 
FOR  NEAR  WHITE  AND  RANDOM  WALK  GENERALIZED  PHASE  MODEL 
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Taking  the  Fourier  transform  of  equation  (4.5a),  we  have  the  power 


spectral  density  of  the  cw  waveform. 


Id 


k  #—  k  cos  b  -  b  sin  b  k  cos  b  -  b  sin  b  - 

2  o  I  o  1 

I  c|  T  - -  exp  {-k  } I  - - - 

o  2  2  o  2  2  2  2 

k+b  »-  k+b  k+b  -J 

—  o  o  1 


2  2 
<T  T  er  T  T 

1  o  2  1  o  o 

where  k  * - +  <r  ,  k  - ■ — ,  and  b  *  w  —  . 

o  4  o  1  4  2 


(4.6a) 


The  power  spectral  density  of  the  generalized  model,  equation 
(4.5b),  which  contains  the  superposition  of  the  "near  white"  stationary 
phase  random  processes  and  the  random  walk  phase  random  process  is 
found  by  taking  its  transform.  The  result  is 


L  r-  exp(-k'  )  T  T 

2  v  J  1-1  r  1-1  1-1 

/  \  - I  exp  {-k - }  [k  cos(  w - 

*—  \  2  2  i  2  i  2 

i“0  k  +  w 


S(w)  -  I  c  | 


)  -w  sin(w  — 
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l  T  T  — w 

1  1  1  -J  1 

-exp{-k  — }  [k  cos(w  — )  -  w  sin(w  — )]  I 

i  2  i  2  2  J  \ 


(4.6b 


where  k'  ,  k  ,  and  T  are  defined  after  equation  (4.5b). 
i-1  i  i 


4.1.2.  The  infinite  pulse  train: 


The  time  average  auto-correlation  function  of  this  envelope  waveform 
is  given  by: 


R(t)=  lim  -  \c  -  >  d(u+t-*T)  c  -  >  d(u-jT)  du 

o  A — >aO  2A  J  M  /\t  —  \l  /\t  ^ 

-A  ^  -o©  -c©- 


where  T  is  the  interpulse  period  and  At  «  1  is  the  pulsewidth.  This 

2 

waveform  has  been  normalized  so  that  the  average  energy  is  I  cl  . 
Setting  i=^-j  and  carrying  out  the  integration,  we  obtain 


i=o© 

2  T  r~ 

R  (t)  =  |  c|  -  y  d(t  -  iT)  (4.7  ) 

o  Ac  *— 

i=-CXD 


For  the  modulated  output  waveform,  we  let  m/2  :  m  even;  denote  the 
nunfoer  of  interpulse  intervals,  T,  over  which  the  white  stationary  term 
is  still  partially  correlated  and  we  allow  decorrelation  to  occur  within 
the  next  interpulse  period.  Then  the  stationary  "near  white” 
decorrelation  interval 


T  T 

o  o 

- <  t  < - 

2  2 

will  be  related  to  the  interpulse  interval,  T,  by 

mT  <T  <(m+2)T 
o 

This  relationship  is  illustrated  in  figure  4.5a.  The  corresponding 

average  auto-correlation  function  for  the  modulated  waveform  is  obtained 

by  multiplying  R  (t)  of  equation  (4.7)  with  R  (t)  of  equation/(  2 . 10 ) 
o  b 

and  is  given  by: 


i-m/2 

2  T  12  2  2 

R(t)  “  |cl  -- —  >  exp { - (  <r+  2  or  — )|  t| }  d(t-lT) 
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I t|  >  To/2 


(4.8) 


Next  we  determine  the  average  power  spectrum  of  the  output  system 
waveform  due  to  the  Infinite  pulse  train.  Here  we  do  not  consider  the 
generalized  phase  model  containing  the  superposition  of  "near  white" 
phase  random  processes  introduced  in  equation  (4.5b),  figure  4.4.  The 
average  power  spectrum  for  such  a  generalized  process  can,  however  be 
found  in  the  same  manner.  Figure  4.5a  shows  the  average 
auto-correlation  function  of  the  product  of  the  infinite  pulse  train  and 
the  oscillator  wav  form.  The  average  power  spectrum  is  given  by: 


2  '  V" 

S(w)  -  |c|  T  I  1+  ^  exp{-ik  T}[exp{jiwT}  +  exp{-jiwT}] 
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o  1  2  1 
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where  a  =  k  T,  b  «  wT ,  k  - - +  <T - , 

o  o  o  2  o  T 
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<r 

1 

- ,  and  a  ■  k  T 
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For  the  special  case,  when  the  near  white  stationary  phase  term 
of  figure  2.1a  becomes  decorrelated  immediately  after  the  first  pulse, 
i.e.  when 

+ 

T  «  T  or  m  -  0 
o 


as  shown  in  figure  4.5b,  the  expression  (given  in  equation  (4.9a))  for 
average  power  spectrum  of  the  infinite  pulse  train  reduces  to 
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(4.9b  ) 


where  a  -  k  T 
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4.1.3.  The  finite  pulse  train  of  N  pulses 


The  time  average  auto-correlation  function  of  the  envelope  waveform  is: 
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where  T  is  the  interpulse  period  and  At  «1  is  the  pulsewidth.  The 
energy  in  the  pulse  train  has  been  normalized  to  I  cl  .  Setting  i“J?-j 
and  carrying  out  the  integration,  we  obtain: 

i-N-1 


2  T  ''T— 

R  (  t  )  ■  I c|  y  (N  -  lil  )  d(t  -  IT) 

o  NA* 


(4.10) 


i— m-i 


Equation  (4.10)  is  plotted  in  figure  4.6. 


FIGURE  4.6  TIME  AVERAGE  AUTO-CORRELAXION  FUNCTION  OF  FINITE  PULSE  TRAIN 
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For  the  modulated  output  waveform  of  the  finite  pulse  train  where 


mT  <  T  <  (m  +  2)T 


the  resulting  average  auto-correlation  function  is 
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Next  we  determine  the  average  energy  spectrum  for  the  finite  pulse 
train.  We  take  the  simplified  case,  m  ■  0.  Taking  the  transform  of  the 
above  average  auto-correlation  function  (equation  (4.11))  we  have  for 
m  *  0,  i.e.  figure  4.5b, 
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(A. 12) 


where  a  -  k  T,  b  «  w  T,  and  k  ■  - 

1  12 


The  average  energy  spectrum  for  more  general  correlation  functions, 
l.e.  for  m  1,  and  where  the  stationary  phase  decorrelation  time 
extends  beyond  the  first  pulse  can  also  be  found  in  a  similar  manner. 
He,  however,  will  not  consider  this  case  here  because  the  expressions 
become  quite  cumbersome  without  adding  any  further  insight  into  the 
method  of  analysis  presented  here. 

Since  the  modulating  envelope  waveforms,  f  (u),  modeled  for  the 

2  o 

oscillator  were  normalized  for  energy  |c|  ,  the  signal  to  noise  ratio 

2 

in  every  case  would  be  |c|  /N  ,  [4.3]  where  N  is  due  to  thermal 

o  o 

additive  noise. 

The  three  equations,  (4.6),  (4.9),  and  (4.12)  provide  the  spectral 
characterization  for  the  following  instability  conditions: 

1.  Equation  (4.6)  -  Unmodulated  CW  oscillator  signal. 

(4.6a)  -  “near  white”  stationary  and  random  walk  [4.4] 

(4.6b)  -  "near  white”  and  intermediate  stationary,  random  walk 
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2.  Equation  (4.9)  -  Infinite  pulse  train  modulated  oscillator  signal 


(4.9a)  -  "near  white”  stationary  and  random  walk,  stationary 

correlation  spanning  more  than  one  interpulse  interval 
(4.9b)  -  stationary  white  and  random  walk,  stationary  correlation 
spanning  less  than  one  interpulse  interval 


1 


r 


1 

\-~ 


1 


% 


3.  Equation  (4.12)  -  Finite  pulse  train  modulated  oscillator  signal, 
"near  white”  and  random  walk,  stationary  spanning 
less  than  one  interpulse  interval. 


4.2 .  Random  Walk  Frequency  and  Frequency  Linear  Drift 

In  this  section,  we  consider  the  second  model ,  where  the  phase 
random  process  is  characterized  by  random  walk  frequency  and  frequency 
linear  drift.  First,  we  treat  the  oscillator  rf  signal  only.  After 
obtaining  expressions  for  the  power  spectral  density,  we  will,  at  the 
end  of  this  chapter,  consider  modulation  by  the  finite  pulse  train. 

The  auto-correlation  function  of  the  oscillator  waveform  has  been 
determined  in  equation  (2.15).  The  power  spectral  density  of  the 
oscillator  rf  signal  having  a  random  walk  frequency  component  can  be 
found  by  evaluating  the  Fourier  transform  of  equation  (2.15). 
Specifically 


« 
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S  (w)  = 
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exp<  - 


-oc 


1 

2  2  2  I 

-  1 1 1  -  jwt  > 


dt 


(4.13) 


pa 


To  the  author's  knowledge,  no  closed  form  or  tabulated  solution  exists 
for  the  above  Integral.  Such  solutions,  however,  do  exist  for  factors 
of  the  auto-correlatlon  function  appearing  In  the  above  Integrand  [4.5]. 
These  can  be  used  successfully  to  evaluate  the  effect  that  random  walk 
frequency  has  on  power  spectral  spread.  The  two  factors  appearing  in 
the  above  integral  give  rise  to  the  random  walk  frequency  along  with  the 
offset  random  frequency.  The  first  factor  (first  term  in  the 
exponential)  is  solely  due  to  the  stationary  independent  increments. 
The  second  factor  (second  tet.n  in  the  exponential)  is  due  to  the  long 
term  frequency  linear  drift,  the  offset  frequency  at  time,  u.  (See 
discussion  on  equation  (C.2)  in  Appendix  C) 

4.2.1.  Frequency  Linear  Drift  Deconvolution 

Equation  (2.15),  for  the  auto-correlation  function 
can  be  written  as  a  product 
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Since  the  product  of  the  auto-correlation  functions  implies  convolution 
of  the  power  spectral  densities  we  have 
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(4.14a) 


S  (w)  -  S  (v)  *  S(w) 
v  2 


where  the  Integral 


Is  the  spectrum  resulting  from  the  offset  frequency  (frequency  linear 
drift)  at  time,  u,  and 


is  the  spectrum  resulting  from  the  random  walk  (stationary  independent 

increments)  frequency  (referenced  to  the  desired  oscillator  frequency  at 

u,  i.e.  t-0).  Note  [4.6]  that  equation  (4.14b)  gives  rise  to  a  power 

spectral  density  with  Gaussian  weighting  over  the  spectral  variable,  w, 

x 

with  mean,  zero,  and  variance,  c^d^  •  This  is  precisely  the  assumption 
we  make  at  time,  u,  for  the  offset  (long  term  linear  drift)  frequency  in 
equation  (2.16)  and  again  in  equation  (C.2b).  In  those  cases,  we  model 
the  random  walk  frequency  waveform  as  having  a  Gaussian  offset 
condition.  Clearly,  from  equation  (4.13)  and  equation  (4.14a),  we  see 
that  the  spectrum  for  the  frequency  linear  drift  component  given  by 
equation  (4.14b)  has  a  dispersive  (broadening)  effect  on  the  spectrum 
for  the  random  walk  frequency  component  given  in  equation  (4.14c).  The 
effect  which  the  offset  Gaussian  distributed  frequency  has  on  the  final 
power  spectral  density  arises  from  the  convolution  given  in 
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equation  (4.14a). 


Time  and  spectral  domain  relationships  are  well  known  [4.5],  [4.6] 
for  the  Gaussian  model  represented  by  equation  (4.14b).  In  the  next 
section  we  will  consider  the  evaluation  of  the  stationary  independent 
Increment  components  of  the  random  walk  frequency,  equation  (4.14c). 

4.2.2.  Series  Solution  for  Related  Airy  Integral 

We  will  focus  attention  on  the  random  walk  frequency  drift  effect 
of  equation  (4.14c).  The  integral  appearing  in  that  equation  is  a 
"related  Airy  function",  properties  and  applications  of  which  have  been 
considered  in  the  literature  [4.5],  [4.7] -[4.11] .  In  order  to  apply  the 
references  in  the  evaluation  of  equation  (4.14c),  we  proceed  as  follows: 
First  we  change  our  variable,  t,  to 


t'  -  -gt  :  t  <  0 
t'  “  gt  :  t  >  0 

Equation  (4.14c)  then  becomes 
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-1  -1 
c©g  OOg 

13  (  13 

exp{ - t  +  jgwt}dt  +  \  exp{ - t  -  jgwt}dt 

3  /  3 

o  o 


] 


(4.15a) 


where  the  complex  character  of  g  Indicates  the  integration  path  which 
will  be  developed  later.  Let  us  write  this  equation  as  follows. 


S(w)  -  TT  g  [  Hi(jgw)  +  Hi(-jgw)  ] 


(4.15b) 


where 
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Hi(x)  A  — —  \  exp{  -  t  /3  +  xt  }  dt  :  x  (real) 

"nr  ; 

-CO 


and 

2  1/3|  2  1/3|  2 

g  »  (2/<T  )  -  j(  2/<r  >  I  exp(j  -  IT  k)  :  k  -  -1,  0,  1 

2  I  2  i  3 


is  In  general,  a  complex  variable.  The  complex  limits  in  equation 
(4.1Sa)  represent  the  radial  path  of  integration,  see  figure  4.7,  and 
Hi(w)  is  the  "related  Airy  function".  For  t  real,  equation  (4.15)  is 
evaluated  by  taking  the  above  integral  over  any  one  of  three  paths  as 
shown  in  figure  4.7.  We  shall  show  that  the  integral  results  in  the 
same  expression  regardless  of  the  path  of  integration,  k  ■  -1,  0,  1, 


taken. 


I 


FIGURE  4.7  INTEGRATION  PATHS  FOR  RELATED  AIRY  INTEGRAL 


From  equations  (6)-(8),  (10)  and  (11)  of  [4.11],  we  can  evaluate 
equation  (4.15)  where  k  *  0,  -1,  and  +1  correspond  to  the  paths  0  to  a©, 
0  to  cud  exp(j2/3  IT),  and  0  to  as  exp(-j2/3  IT)  respectively.  Let  us 
rewrite  equation  (4.15)  with  the  following  substitution; 


z  -  jgw  =  | w | exp ( j  TT/2)  |g|exp(j  2kIT/3) 

=  I z |exp( j  [IT/2  +  2krT/3]) 


:  k  *  0,  -1,  1 


which  preserves  the  integration  paths  of  figure  4.7.  Then  equation  (4.15) 


is  evaluated  using  equations  (10)  and  (11)  of  [4.11],  where  1,1,  and 

1  2 


I  in  that  reference  are  designated  by  k  =  0,  -1,  and  1  respectively.  We 
3 

then  have 


o o 
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Clearly  for  odd  values  of  r  the  above  equation  is  zero.  Summing  over 
even  values  of  r,  multiplying  g  by  the  first  exponential  in  the  above 
equation  and  substituting  the  above  expression  for  z  we  have 
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which  is  invariant  with  k  *0,  -1,  and  1.  Furthermore 


Ue  may  therefore  take  the  integral  over  the  path  0  to  oo,  i.e.  k  witl 


1/3 


in  order  to  determine  the  spectrum,  S(w). 


4. 2. 2.1.  Limits  of  the  Series  Solution 


The  series,  given  in  equation  (4.16),  is  useful  for  small  values  of 
w.  When  using  equation  (4.16)  to  compute  S(w),  it  must  be  kept  in  mind 
that  the  maximum  term  in  the  series  will  be 


1/3 

1/3  3  p 

-2/3  (3  p) 

max{s }  -(2)3  -  :  p  -  |z  I 

1/3 

(3  p)l 


Solving  this  expression  for  p,  we  find  the  index,  r,  which  gives  the 
largest  term  in  the  series.  This  occurs  at 


r 


1/3 
3 

- p 

(max)  2  (max) 


It  is  therefore  necessary  that  equation  (4.15)  be  used  only  for 


w  <  p 

(max) 


z(max) 


in  order  to  be  within  the  precision  limits  of  the  computational  device. 

4.2.3.  Asymptotic  Expansions  for  Related  Airy  Integrals 
For  cases  where 


w  >  p 


2 


1/3 


asymptotic  expansions  for  equation  (4.16)  can  provide  results  well 
within  most  precision  requirements.  From  Poincare  [4.12],  Scorer  [4.7], 
and  Olver  [4.13]  the  following  asymptotic  expansion  is  obtained  for 
sufficiently  large  negative  real  x 


1 

Hi(x)  "  - 

rrx 


an 


V 


1  + 


(3s  +  2)! 


3  s 

x  ‘ -  s  !  (  3x  ) 

s-0 


x  — ^  -ao  (4.17) 


where  the  symbol  ~  denotes  "asymptotic  to”.  The  error  of  the  above 


expansion  continues  to  decrease  so  long  as 


3 

(3s+2)(3s+l)(3s)/s  <  3 | x | 

3/2 

s  <  < 1/3) lx! 

for  large  negative  x,  say  |x|>5.  In  the  complex  plane,  the  asymptotic 
expansion  -  with  x  (real)  replaced  by  z  (complex)  -  holds  [4.13]  for  the 
analytic  continuation  Hi(z)  of  Hi(x),  provided  that 

|arg(-z) |  <  IT/3  -  d, 

(it  holds  for  2  IT/ 3  +  d  _<  arg(z)  _<  4  IT/ 3  -  d:  it  does  not  hold  for 
-2  IT/3  +  d  _<  arg(z)  <  2  IT/3  -  d:  d  ■  an  arbitrary  positive  constant.) 

We  will  next  modify  the  expression  in  equation  (4.17)  so  that  it 
will  be  analytically  contlnous  along  the  jw  axis.  We  use  a  method 
developed  by  Olver  [4.13]  which  assures  analytical  continuity  along  the 
imaginary  axis  through  a  rotation  of  the  coordinate  system. 

4. 2. 3.1.  Analytic  Continuation  at  z  *  £jw 
When  evaluating  equation  (4.17)  for 

2  *  ±jgw  "  gw  exp(-jB)  :B  -  (IT/2) 

we  apply  Olver  [4.13],  theorems  3.2  and  3.3  to  the  integral,  Hi(jjgw), 
i.e.  chapter  11,  equation  (12.12)  [4.13],  and  use  the  principal  value  of 


Watson's  Lemma,  [4.13],  states  for  certain  convergence  conditions,  that 
the  Integral 


on 

-zt 

I(z)  -  \  e  q(t )  dt 


has  an  asymptotic  expansion 


oo 


■<e —  | - >  s  +  A  i 

>  <  -  >  - 

Z —  |  u  (s  + 


I(z)  ~ 


A)/u 


s-0 


(8  +  A)/u 

as  z  — >  oo  in  the  sector  |arg(z)|  <  IT/2  where  z  has  its 

principal  value  and  q(t)  has  the  Mclaurin  expansion  as  t  — >  CH- 


q(t) 


o© 

(s  +A  -  u)/u 

a  t 
s 

s«0 


c.f.  theorem  3.2,  page  113,  [4.13].  For  the  case  of  interest  the  para¬ 
meters,  A*u*l*  From  theorem  3.3,  page  114,  [4.13],  and  the  discussion 
on  page  431  [4.13]  we  may  find  the  analytic  continuation  of  I(z)  in  the 
sector,  -  2  IT/3  <  arg(z)  <  2  n/3  by  a  rotation  of  the  path  of 
integration  by  B  n/6.  Applying  theorem  3.2,  [4.13],  with  t  exp(jB) 
and  z  exp(-jB)  playing  the  roles  of  t  and  z  respectively,  we  have  for 
the  Mclaurin  expansion,  equation  (3  .02),  page  113,  [4.13], 


(4.18a) 


jB  1  3  j3B 

q(te  )  ”  exp( -  t  e  ) 

3 


2  j3B  jB  3  40  j6B  jB  6 

1 - e  (te)+  —  e  (t  e  )  + 

3!  6! 


We  see  that  the  constants  in  the  Mclaurin's  expansion,  page  113, 


equation  (3  -02),  [4.13],  become 


2  j3B  40  j6B  2240  j9B 

a  =1,  a  =  —  e  ,  a  ”  —  e  ,  a  ■  — —  e  ,  etc. 
o  3  3!  6  6!  9  9! 


The  asymptotic  expansion,  equation  (3  .03),  page  113,  [4.13],  becomes 


1  - »  2 

I(z)  -  -  -  (4)  — 

-jB  3! 


2  j3B  1 


4  “j4B 


r40  j6B  1 

(7)  —  e  - 

6!  7  -J7B 


1  - »  444U 

-  -  (10)  - 

7  -J7B  9! 


2240  j9B  1 


+  ....  (4. 18b ) 


10  -jlOB 


In  order  to  evaluate  the  integral,  I(z)  at  z  =  £jw  -  w  exp(£jIT/2), 
we  rotate  the  paths  of  integration  by  setting  B  -  ^(IT/6)  respectively 
for  the  £  operator,  i.e. 


l(±j«)  “  I(jv;  B  -  n/6)  or  I(-jw;  B  *  -IT/ 6) 


We  obtain 
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i  ±jn/6  2  jjrr/6  40  jjrx/6  2240  jjrr/6 

I(+jw)  =  +  —  e  +  j  —  e  T -  e  -  j -  e  ±  •  •  • 

4  7  10 

jw  w  jw  w 


4  7  10 

I(jtJw)  =  (Jl/jw  +2/w  T  40/ jw  -2240/w  +  ...  )  exp(j-jIT/6) 


(4.18c) 


Combining  the  ^  components  as  the  superposition  of  the  two  expressions 
for  <z,B>  =•  <+jw,rT/6>,  <-jw,-lT/6>,  we  have 


I(jw)  +  I(-jw)  =  (  1/w  -40/w  +  ...)  +  \J"3  (  2/w  -2240/w  +  . ..)  (4.18d) 


Substitution  of  the  parameters  from  equation  (4.15)  results  in  the  power 
spectral  density  asymptotic  expansion. 


S(w)  “  g 


w  (gw) 


1  40  2\J3  2240\j3 


4  10 

(gw)  (gw) 


roc 
1  «$- — 

s<w)  -  ;  <J  i  -  — 


s  (6s+5)! 


6  c —  s  68 

3(gw)  s-0  ( 2s-*-l) !  9  (gw) 


{(ED 

g  w  s»0 


s  (6s+2) ! 

I  — — 

s  6s 
(2s)!  9  (gw) 


(4.19) 


where 


2  1/3 

g  -  (2  /  <r  ) 

2 


-  102  - 


« 

and  w  is  sufficiently  large.  This  is  the  asymptotic  expression  for  the 

power  spectral  density  of  a  cw  rf  signal  having  a  random  walk  frequency 

2 

m  instability  with  drift  variance  O'  • 

V  2 


N 

|T 

■ 

L 

% 


4.2.4.  Effects  on  Modulating  Waveforms 

Next  we  consider  the  two  modulating  waveforms,  the  cw  and  the 
finite  pulse  train. 

4. 2. 4.1.  CW  Waveform 

When  the  modulating  waveform  is  a  cw,  with  amplitude  c,  we  have  for 

2 

power  spectral  density,  the  product  of  equation  (4.19)  and  |c|  ,  a 
constant. 


2 

S  (w)  -  |c|  S(w)  (4.20) 

cw 

4. 2. 4. 2.  Finite  Pulse  Train 

The  time  average  auto-correlation  function,  Ro(t),  of  the  finite 
pulse  train  envelope  was  previously  given  on  equation  (4.10).  The 
auto-correlation  function  of  the  rf  oscillator  signal  corrupted  by  the 
random  walk  frequency  instability  onJy,  i.e.  neglecting  any  offset 
frequency  Instability,  is  found  in  equation  (2.15) 


R  (t) 
b 


2  3 

<r  Itl 
2 


6 


(4.21) 
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Hence,  the  average  auto-correlation  function  of  the  modulated  output 
waveform  of  the  finite  pulse  train  is 


i-(N-l) 


2  T 

R(t)  -  R  (t)R  (t)  -  |c |  - 

o  b  N/\t 


>  d(t-iT) 


i— (H-l) 


(4.22) 


Ue  determine  the  average  energy  spectrum  of  the  random  walk 
frequency  rf  signal  modulated  by  the  finite  pulse  train  by  taking  the 
Fourier  transform  of  equation  (4.22).  The  result  is 


(4.23) 


A  closed  form  or  tabulated  solution  for  equation  (4.23)  was  not 
attempted  due  to  the  analytical  difficulties  encountered  from  that 
expression.  A  numerical  evaluation  will  be  accomplished  in  following 
chapter. 
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V  NUMERICAL  EXAMPLES  AND  APPRO’  IMATIONS  FOR  POWER  SPECTRAL  DENSITY 

In  the  previous  chapter,  we  have  considered  the  spectral  spread  due 
to  two  groups  of  phase  random  processes.  Throughout  this  chapter,  we 
present  numerical  results  obtained  for  spectral  spread.  Since  the 
spectral  representations  are  symmetrical,  we  will  present  the  one  sided 
response  only.  For  some  cases,  approximations  were  required  and 
approximate  results  are  presented.  As  we  have  done  thus  far,  the  two 
models  will  be  treated  separately. 

5.1 .  White  Phase  and  Random  Walk  Phase  -  Oscillator  RF  Signal  with 
Modulation  Waveforms 


5.1.1.  CW  Waveform 

Suppose  a  CW  oscillator  has  a  white  phase  component  with  standard 
deviation 

-10 


<r  -  10  radians 


o 

over  an  averaging  time  of 


seconds 


and  a  random  walk  phase  component  with  standard  deviation 


-9 

<r  -  io 

l 


Then  the  auto-correlation  function  of  the  output  waveform  is  shown  in 
figure  5.1a  and  the  resulting  one  sided  power  spectral  density  obtained 
from  equation  (4.6a)  is  shown  in  figure  5.1b. 
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When  O’  ,  <r  ,  and  T  are  small,  i.e.  «  1,  then  the  power  spectral 
o  1  o 

density  for  the  cw  oscillator,  equation  (4.6a),  can  be  approximated  over 
the  following  frequency  Intervals: 


1  1 

-  S(b)  -  -  :  b  «  It  <  k  «  IT  (5.1a) 

2  k  1  o 

I  c  |  T  1 


Definitions  for  the  parameters  k  ,  k  (k  <  k  ),  and  b,  will  be  found 

olio 

immediatly  following  equation  (5.2e). 


In  Rb(  t) 


(a)  Example  -  Auto-Correlation  Function  for 
White  Phase  &  Random  Walk  Phase 


FIGURE  5.1  EXAMPLE  -  CW  WAVEFORM 


al  Density  of  CW  Oscillator 


CW  WAVEFORM 


Equation  (5.1a)  Is  obtained  through  the  cancellation  of  the  first  two 
terms  of  equation  (4.6a),  with  the  approximation  of  its  last  term 


f r 


► 


remaining.  Equating  b  to  k  ,  the  last  term  of  equation  (4.6a)  becomes 

1 

1 

- -  :  b  -  k  <  k  «  IT  (5.1b) 

2k  1  o 

1 


When  b  is  in  the  range  between  k  and  k  ,  still  very  small,  the  trigono- 

1  o 

metric  functions  in  the  numerators  of  equation  (4.6a)  are  approximated 
by  one  and  zero  for  the  cosine  and  sine  respectively.  The  last  becomes 


k 

1  :  k  <  b  <  k 

-  1  o  (5.1c) 

2  2 

k+b  :  k  <  b  «  n/  2 

1  o 


Equation  (5. Id)  is  obtained  as  follows:  With  k  and  k  very  small 

1  o 

the  bracketed  part  of  equation  (4.6a)  is  written 


(1-k  ) 
o 


k  cos  b  -  b  sin  b  -  k  cos  b  +  b  sin 
o  1 


b  ] 


— . 

_ 


setting  1-k  ~  1,  cancelling  terms  and  factoring  the  result,  we  have 


k  -  (k  -  k  )  cos  b 
o  o  1 

-  :  k  «  b  <  n/2  (5. Id) 

2  o 

b 


When  b  is  equal  to  IT/2  it  dominates  all  other  terms  in  equation  (4.6a). 
The  cosine  terms  are  zero,  the  sine  terms  are  one,  and  we  have 
1 

— —  (  k  -  b  +  b  )  which  leaves  us  with 
2  o 
b 
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i 


o 


:  b  -  IT/ 2 


(5.1e) 


k 

2 

b 


In  the  above  expressions  we  have  made  the  following  substitutions. 


cr  T 

1  o  2 

k  * - +  cr  , 

o  4  o 


Note  that  other  notations  [5.1] 
parameters 


2 

cr  T  T 


[5.2] 


1  o  o 

- —  ■  <  k  ,  and  b  *  w  — . 
4  o  2 

can  be  related  to  the  above 


T  h  T 

o  2  o 

k  ■  h  -  +  - -■ ,  k  «■  h  - 

o  o  4  2  1  o  4 

(2  IT) 


Figure  5.2a  shows  the  relationship  of  the  parameters  in  equation 
(5.1)  to  the  auto-correlation  function  for  the  cw  oscillator  signal. 
Figure  5.2b  illustrates  the  plotting  of  the  power  spectral  density 
curve . 


5.1.2.  Infinite  Pulse  Train 

We  can  approximate  equation  (4.9b)  and  obtain  the  powei  spectral 
density  when  the  stationary  component  is  white.  Approximations  to  the 
infinite  pulse  train  example  are  made  as  follows: 


111  - 


In  Rf  =  In  |c|2+lnRb(t) 

I  I  V* 


T 


N 


v 


FIGURE  5-2  DEPENDENCE  OF  AUTO- CORRELATION  FUNCTION  AND  POWER  SPECTRAL 
DENSITY  ON  PHASE  STATISTICS  -  CW  WAVEFORM 
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Given,  the  average  power  spectrum  for  the  infinite  pulse  train, 
equation  (4.9b),  and  denoting  its  left  side  as  a  normalized  term  we  have 


S(w)  2  cos  b  -  exp(-a) 

S  (w)  - - *  [  1  +  exp(-<r  )  - 

N  2  o  cosh  a  -  cos  b 

I  c  |  T 


where  a  «  1  and  <T  «  1.  For  the  following  conditions:  S  (w)  is 

o  N 


2  1  -  (  1  -  a  )  2  2  2 

b  -  0:  “  [  1  +  (1  -  ff) - ]  I  1  +  (  1-0 - 2  — 

o  2  o  a  a 

1  +  a  /2  -  1 


2  cos  a  -  exp(-a) 

b  «  a:  2  [  1  +  (1  -  <T  ) -  ] 

o  cosh  a  -  cos  a 


2  2 
2  l-a/2-(  l-a+a/2) 

[  1  +  <1  -  O - ] 

o  2  2 

1  +  a  /2  -  1  +  a  /2 


"  1  +  (  1  - 


2  a  -  a 

O'  ) - 

o  2 


"  1  +  (  1/a 


-  °J‘  - 


1  +  <r  )  ~  1/a 


a  «  b  <  tt: 

2  cos  b  -  exp (-a) 

2  1+  exp(  -  O'  ) - 

o  cosh  a  -  cos  b 


2  cos  b  -  (1  -  a) 

1  +  (  1  -  ff  ) - 

o  1  -  cos  b 


2  a 

1  +  (  1  -  cr  )  ( - 1  ) 

o  1  -  cos  b 


1  + - 1 

1  -  cos  b 


a  cr 

o  2 

1  -  cos  b  o 


2  a 

<r  + - , 

o  1  -  cos  b 


The  approximations  are  thus  summarized  as  follows: 


1  2 

-  S(b)  -  — 

2  a 

|c  |  T 


b  -  0 


:  b  ■  a 


2  a 

c r  + - 

o  1  -  cos  b 


:  a  «  b  <  IT  (5.2) 


where  a  - 


k  T,  k  - 
1  1 


and  b  -  wT . 


Equation  (5.2)  is  sketched  for  0  <  b  £  IT  in  figure  5.3.  Note  that  the 
peak  of  the  response  is  limited  by  the  random  walk  component  through  the 
parameter,  a.  The  noise  floor  (null  depth)  is  determined  by  the  white 
phase  component,  the  random  walk  phase  component,  and  frequency. 
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Oscillator  & _ 

julse  Train  Spectrum  Dominated 


FIGURE  5.3  ONE  SIDED  AVERAGE  ENERGY  SPECTRUM  -  INFINITE  PULSE  TRAIN 


5.1.3.  Finite  Pulse  Train 

Suppose  a  train  of  N*32  pulses  with  Interpulse  period  T-.001 
modulates  an  oscillator.  Let  the  oscillator  have  a  white  phase 
component  with  standard  deviation 

-8 

<T  -  10 
o 


over  an  averaging  time  less  than  one  interpulse  period,  T,  and  a  random 
walk  phase  component  with  standard  deviation 


-7 


O' 

1 

-  10 

radians/  \jsec 

The 

pulse 

train 

is  shown  in 

figure 

5.4a,  the  auto-correlation 

function 

of  the 

output 

waveform  is 

shown 

in  figure  5.4b,  and  the 

resulting  average  energy  spectrum  obtained  from  equation  (4.12)  is  shown 
in  figure  5.4c.  Table  5.1  gives  the  peaks  and  nulls  of  the  average 
energy  spectrum  from  b”0  to  b^TT.  In  plotting  the  spectrum  of  the 
finite  pulse  train,  we  observe  that  a  number  of  approximations  to 
equation  (4.12)  can  simplify  the  computations.  These  simplifications 
allow  us  to  determine  the  peaks  of  the  spectrum  sidelobes  and  the  null 
depths,  and  are  written  below: 


2 

c r  T 
1 

For  0  <  d  ,  <T  «  1,  a  - - ,  b  *  wT  we  have  for  the 


o  1  2 

1 

Mainlobe  - S(b)  =  N  :  b  -  n  IT 

Peak  2  :  0  <  n  even 


I  c  |  T 


Lowest  1  2 

Peak  - S(b)  =  cr  +  a 


2  o 

I  c  |  T 


b  =  n  IT 
n  odd 


2 

Sidelobe  - -  :  b  *  n  IT  /  N 

Peaks  N( 1  -cos  b) 

:  0  <  n  <  N;  n  odd 


Null 

Depth 


2  2a 

*=  <T  + -  :  b  =  n  TT  /  N  (5.3) 

o  I  -cos  b 

:  2  <  n  <  N;  n  even 
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la)  n  =  32  Finite  Pulse  Train 


CORRELATION  FUNCTION 


lb)  CORRELATION  FUNCTION  OF  PRODUCT  OF  FINITE  PULSE  TRAIN, 
NEAR  WHITE  PHASE  (m-O).  and  RANDOM  WALK  PHASE 


FIGURE  5.4  PROPERTIES  OF  MODULATED  32  PULSE  TRAIN 


m  I  ttt 
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29  ?<  ?5  ?7  ?«  79  30  3]  V, 


32  b/rr 


(c)  One  Sided  Average  Energy  Spectrum  -  White  &  Random  Walk  Phase 
FIGURE  5.4  PROPERTIES  OF  MODULATED  32  PULSE  TRAIN 


The  spectrum  sidelobe  peaks  and  null  depths  (upper  and  lower 


spectrum  limits)  are  plotted  in  figure  5.5  for  various  parameters  as 
Indicated.  These  plots  are  for  uniform  weighted  pulse  trains.  With 
some  effort  it  is  a  straight  forward  procedure  to  obtain  the  spectrum 
for  a  weighted  train  of  pulses.  Mote  from  equation  (5.3)  and  figure  5*5 
that  for  the  finite  pulse  train  the  peak  response  is  determined  by  the 

number  of  pulses  and  frequency,  w;  and  the  null  depths  are  determined  by 

2.  k 

the  white  phase  component  (Jo,  the  random  walk  phase  component,  01,  and 

the  frequency,  w  for  all  but  large  values  of  the  phase  variance.  When 

the  variance  is  large,  then  the  spectral  spread  deteriorates  beyond  that 

of  the  sidelobe  and  mainlobe  peaks  as  shown  in  figure  5*5b. 


5.2 .  Random  Walk  Frequency  -  Oscillator  RF  Signal  with 
Modulating  Waveforms 

To  evaluate  the  effects  of  spectral  spreading  of  an  oscillator  when 
random  walk  frequency  instability  is  present,  we  use  the  exponential 
expansion,  equation  (4.16),  and,  the  asymptotic  expansion, 
equation  (4.19),  in  normalized  form.  We  perform  the  normalization  by 
setting  the  argument,  gw*ci,  where 


g  * 


1/3 


We  obtain 

GO  I  1/3  1 2r 

1  -2/3^—  J  3  ci)  r 

—  s  (w)  -  2  (3)  >  -  (-1) 

g  v  ^ —  ( 2r ) ! 

r-0 


2 


-  [r-1] ) ! 
3 


(5.4a) 


for  the  exponential  expansion  and 


1  1 
—  s  (w)  -  — 
g  v  ci 


3 

(6s+2) !  (ci)  -  (6s+4)(6s+5)] 

8  6s+6 

(2s)J  9  (ci) 


(5.4b) 


for  the  asymptotic  expansion.  Figure  5.6  illustrates  the  behavior  of 
the  oscillator  power  spectral  density  with  normalized  frequency,  i. 


For  the  finite  pulse  train,  the  average  energy  spectrum  is  found 


through  a  numerical  evaluation  of  the  Fourier  Transform  and  is  shown  in 
figure  5.7. 
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POWER  SPECTRAL  DENSITY  OF  OSCILLATOR 
AND  RANDOM  WALK  FREQUENCY  INSTABILITY 


FIGURE  5.6  POWER  SPECTRAL  DENSITY  -  RANDOM  WALK  FREQUENCY  INSTABILITY 
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FIGURE  5.7  AVERAGE  ENERGY  SPECTRUM  -  RANDOM  WALK  FREQUENCY  & 
32  FINITE  PULSE  TRAIN 


We  normalize  equation  (4.23)  by  setting 

3  2  / 

vT  =  IT  k/N  and  K  -  T  cr  /  2 


Hence 


1  Y"  N  -  i  1  It 

s(k)  -  1  +  2^  - -  exp( - Ki)  cos(n—  i)  (5.5) 


2 

I  c  |  T 


For  K>0,  the  average  energy  spectrum  density  is  that  of  a  perfectly 
stable  pulse  train.  Table  5.2  shows  the  percent  deterioration  of  the 
average  energy  spectrum  for  various  values  of  the  parameter,  K. 

The  parameter  used  by  a  number  of  investigators  [5.1],  [5.2],  [5.3] 


h  «  <r 
-2  2 

2 

corresponds  to  the  frequency  drift  variance,  <r  ,  adopted  here. 

2 


TABLE  5.2 

A  %  DETERIORATION  OF  POWER  SPECTRAL  DENSITY  WITH  PARAMETER,  K 


MAINLOBE 

NEAR 

SIDELOBES 

FAR 

SIDELOBES 

>  -10 

>  -10  “5 

m 

• 

O 

7 

A 

-0.054 

-0.  26 

-0.29 

-27.67 

-54.2 

-49.95 

-59.18 

ROUND  OFF 
ERRORS 

ROUND  OFF 
ERRORS 
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VI  SUMMARY  AND  CONCLUSIONS 


6.1 .  Summary 

We  have  developed  a  procedure  for  finding  the  degradation  in 
spectral  resolution  caused  by  the  phase  instability  of  the  system 
coherent  oscillator.  Phase  instability  was  modeled  as  a  superposition 
of  a  variety  of  random  processes  including  a  "near  white”  Gaussian 
process,  a  non-s tat ionary  phase  and  frequency  random  walk  (Wiener) 
process,  instabilities  at  some  reference  time,  and  long  term  drift 
effects.  The  procedure  involved  determining  the  covariance  matrix  of 
the  phase  process,  obtaining  the  auto-correlation  function  and 
transforming  to  the  power  spectral  density.  Modulated  signals,  i.e., 
signals  obtained  by  modulating  the  oscillator  output  signal  by  some 
deterministic  waveform  were  also  included.  Three  specific  waveforms 
were  treated:  cw,  infinite  pulse  train  and  finite  pulse  train.  Closed 
form  expressions  were  derived.  Approximations  were  used  for  low  level 
phase  instability.  In  our  work,  the  pulse  width  of  the  single  pulse  is 
assumed  to  be  much  smaller  than  the  interpulse  period  so  that  its  effect 
may  be  neglected.  The  results  obtained  here  may  be  applied  to 
communication  and  radar  systems. 

In  addition,  stationarity  and  ergodicity  properties  of  the 
oscillator  rf  signal  containing  the  instability  properties  were  examined 
and  results  were  presented  in  terms  of  a  number  of  theorems. 


6.2.  Conclusions 


The  following  conclusions  are  drawn  from  our  work. 

1)  The  auto-correlation  function  of  the  oscillator  rf  signal  is 
determined  by  evaluating  the  characteristic  function  of  the 
phase  random  process  at  appropriate  values  of  the  argument. 

2)  For  random  walk  frequency,  the  power  spectral  density  of  the 
rf  signal  is  determined  through  evaluating  a  "related  Airy 
integral"  in  the  appropriate  region  of  the  complex  ploie.  For 
typical  values  of  the  drift  frequency  variance,  i.e.  very 
small,  asymptotic  expansions  are  useful  for  even  nominal 
values  of  frequency. 

3)  The  rf  signal  as  effected  by  random  walk  frequency  instability 
is  not  ergodlc. 

4)  The  expressions  for  power  spectral  density  derived  herein 
establish  the  lower  end  of  the  system  dynamic  range. 

o  The  variance  of  the  white  phase  component  establishes  a 
lower  limit  for  the  null  depths  of  the  power  sp-t  :tral 
density  of  the  infinite  pulse  trains  considered.  For  cw 
oscillators,  this  variance  in  conjunction  with  the  random 
walk  variance  establishes  the  spurious  spectral  level  for 
large  frequency  values,  w. 
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o  For  small  values  of  frequency,  w,  the  variance  of  the 
random  walk  phase  component  establishes  the  peak  value 
and  the  roll-off  rate  for  the  power  spectral  density  of 
the  infinite  pulse  train  and  the  cw  oscillator.  This 
variance  in  conjunction  with  the  variance  of  the  white 
phase  component  established  the  upper  limit  for  the 
spurious  spectral  level  of  the  infinite  and  finite  pulse 
train. 

o  For  small  levels  of  instability,  the  number  of  pulses  in 
the  finite  pulse  train  determines  the  peak  main  response 
and  peak  side  lobe  response  of  the  power  spectral 
density.  When  the  instabilities  become  very  large,  i.e. 
the  variances  increase  without  bound,  then  the  side  lobe 
level  increases  with  a  complete  deterioration  of  the 
pattern  including  the  main  lobe. 


6.3.  Future  Work 

The  expressions  generated  herein  have  treated  white  and  random  walk 
phase  separately  from  random  walk  frequency.  Hence,  the  assessment  of 
oscillators  when  subjected  to  internal  or  external  white  noise  sources 
on  one  end  or  environmental  conditions  on  the  other  end  can  be  achieved. 
Expressions  representing  the  composite  case;  white  phase,  random  walk 
phase,  and  random  walk  frequency  is  a  worthwhile  extension  of  these 
results  and  deserves  further  work. 
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Computer  simulations  of  these  phase  random  processes  would  verify 
the  results  and  lead  to  an  experimental  approach  for  evaluating 
oscillators  in  the  laboratory. 

Modeling  of  flicker  noise  should  be  considered  from  the  point  of 
view  of  physical  phenomenon  as  was  done  with  the  phase  instability 
models  introduced  here. 

An  investigation  into  the  synchronization  and  timing  of  multiple 
interconnecting  systems  should  be  performed  using  the  oscillator 
instability  models  developed  in  this  dissertation. 

The  work  performed  in  this  dissertation  can,  with  some 
modifications,  be  applied  to  the  modeling  of  random  effects  in  antenna 
apertures.  The  analogies  between  time -frequency  and  aperture-pattern 
can  then  be  used  to  evaluate  the  limitations  which  random  effects  have 
on  the  performance  of  antenna  systems. 
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APPENDIX  A 


A.O.  Characteristic  Functions  for  Determining  Moments 


of  the  RF  Signal,  b(u). 


The  characteristic  function  can  provide  any  order  moment  [A.l]  of 

the  signal,  b(u),  for  any  distribution  of  the  phase  random  vector,  z- 

For  example,  the  second  and  fourth  order  moment  required  often  in  the 

text  will  now  be  evaluated  for  the  case  when  the  random  phase  vector,  £, 

is  joint  Gaussian  distributed.  From  Davenport  and  Root  [A. 2],  pp  54  and 

152,  equation  (8-44),  the  Gaussian  joint  characteristic  function  is  F  (p) 

£ 

N 

F  (jp  ,  •  ••  ,  Jp  )  ^  E{exp(  j  \  p  z  )]  -  E[exp(j£  2  )] 
z  1  N  —  n  n 


exp (-1/2  £  M£) 


(A.l) 


where  £  is  a  joint  Gaussian  random  vector 


z(u+8+t) 

z(u+t) 

z(u+s) 


with  covariance  matrix 


M  -  E[  £  z  J  -  [  m  ] 
z  ij 


s  1 ,  1  1,  2,  3,  4,  .. 


To  determine  the  fourth  order  moment. 
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* 


R  (t)  -  E[  b(u+s+t)b  (u+t)b  (u+s)b(u)  ]  (A. 2) 


we  proceed  to  evaluate  equation  (A.l)  for  N=4  with 


P  =  1>  P  ”  -1,  p  =  -1,  p  *  1. 
12  3  4 


In  vector  notation  this  is 


Equation  (A. 2)  is  evaluated  using  equation  (A.l)  with  appropriate 
substitutions  for  p,  N,  M,  and  the  random  variable  z.  We  have 


R  (t)  -  F  (j,-j,-j,  j)  &  E[exp( jz  -  jz  -jz  +  jz  )] 
z  12  3  4 

(A. 3a) 


■  exp(  -l/2(m  +  m  +m  +m  }  +  (m  +m  ■+•  m  +  m  )  -  (m  +  m  )) 

11  22  33  44  12  13  24  34  14  23 


In  a  similar  manner  the  second  order  moment  is 


R(s)  -  F  (j,-j)  A  E[exp( jz  -  jz  )]  -  exp(-l/2(m  +  m  )  +  m  )  (A. 3b) 
z  12  11  22  12 


other  moments  of  b(u)  -  exp(Jz(u))  can  generally  be  determined  by  the 

* 

method  just  presented  by  assigning  p-1  for  b(.)  and  p»-l  for  b  (.). 
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APPENDIX  B 


DERIVATION  OF  THE  COVARIANCE  MATRIX  FOR  COMBINED 
WHITE  PHASE  AND  RANDOM  WALK  PHASE  INSTABILITY,  x  +  y 


The  covariance  matrix,  M,  of  the  zero  mean  random  vector,  z 
is  defined  as  the  expectation 


Since  x  and  y  are  independent  we  have 


For  the  Wiener  component,  x,  the  covariance  matrix  is 


id 


3 


Var{x(u+t)}  Cov{x(u+t)  x  (u)} 


Cov{x(u)  x  (u+t)}  Var{x(u)} 


(B.3) 


Substituting  equation  (2.4)  of  the  text  into  the  above  we  have 


(c  +  1 1 1 )  c r 
1  1 


c  <x 
1  1 


2 

c  <r 
1  1 


c  <r 
1  1 


(B.4) 


In  a  similar  manner  we  can  write  the  covariance  matrix  for  the 
"near  white"  component,  y(t). 


Var{y(u+t)}  Cov{y(u+t)  y  (u)} 


Cov{y(u)  y  (u+t)}  Var(y(u)} 


(B.5) 


Substitution  of  equation  (2.3)  of  the  text  into  the  above  results  in 


c  O’ 


2  2 

(c - 1 1 1 )  O’ 

o  To  o 


2  2 

(c - 1 1 1 )  O' 

o  To  o 


c  O' 


:  l t |  <  —  (B.6a) 

2 
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APPENDIX  c 


tr- 


DERIVATION  OF  COVARIANCE  MATRIX  FOR  PHASE  INSTABILITY,  v,  DUE  TO 
RANDOM  WALK  FREQUENCY  AND  LONG  TERM  FREQUENCY  LINEAR  DRIFT 

C.O.  Introduction 

In  this  appendix  we  will  derive  the  covariance  matrix  of  the  phase 
random  process,  v(u),  resulting  from  the  integral  of  a  random  walk 
frequency  process,  v(u),  i.e., 

u+t 

v(u+t)  -  ^  v(y)  dy  +  v(u).  (C.l) 

u 

For  verification  purposes,  we  will  use  three  different  approaches 
to  arrive  at  the  same  result.  We  assume  that  the  random  walk  frequency 
process,  v(u)>  is  a  zero  mean  Wiener  process  with  covariance  function 


.  .*T 
M.  -  E{v  v_  } 
v 


2 

(c  +  1 1 1  )<r 
2  2 


2 

c  o’ 
2  2 


2 

c  <r 
2  2 

2 

c  <r 
2  2 


(C.2a) 


where 


v(u+t) 

v(u) 


m.  *  E{ v}  ■  0^, 
v 
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o*  is  the  random  walk  frequency  drift  variance  and  c  <S  is  the 
2  .  2  2 
frequency  offset  variance,  var{v(u)},  due  to  (past)  long  term  linear 

drift  [C.l],  [C.2]  at  reference  time,  u.  We  assume  that  the  phase 

random  process,  v(u),  also  has  a  zero  mean  Gaussian  random  phase, 

independent  of  v(u),  at  reference  time,  u,  with  an  arbitrary 
2 

variance,  c  <r  .  These  conditions  may  be  stated  as  follows: 

3  2 


.*  2 

Var{v(u)}  *  E{v(u)  v  (u)}  ■  c  o' 

2  2 

*  2 

Var{v(u)}  =  E{v(u)  v  (u)}  *  c  <T 

3  2 


(C.2b) 


(C.2c) 


C.l.  Approximation  Function  Approach 

We  shall  first  use  a  heuristic  approach  based  on  a  direct 
application  of  the  approxinating  function  [C.3]  for  determining  the 
variance ,  Var {v(u+t ) } . 

Suppose  we  assume  the  frequency  takes  on  discrete  random  jumps, 

v(u  ):  u  <  u  <  u  <  «.  <  u  <  «.  <  u  according  to  the  random  walk 
i  o  1  2  i  n 

process.  The  frequency  may  also  be  written  as  a  summation  of  independent 


terms 


V  “V  +(v  -v  )  +  (v  -  v  )  +  . . .  (v-v  ) 
i  o  1  o  2  1  i  i-1 


k-i 


k-1 


where  the  parenthesized  quantities  are  the  k  th  Independent  Increments 

•  •  • 

and  a  notational  change,  v  -v(u  )  was  made,  v  is  the  offset  random 

11  o 

frequency . 

The  phase  random  process,  v(u),  is  the  Integral  of  the  above 
frequency  process,  the  discrete  form  of  which  may  be  written: 


n 

v  -  v  +  >  v  (u  -  u  ) 
n  o^—l  1+1  1 

i-o 


■  v+  v(u-u)+v(u-u)+v(u-u)+.«v(u  ~u) 

oolo  121  232  n  n+1  n 


i-n 


i-o 


(u  -  u  ) 
1+1  1 


(C.4) 


where  v  is  an  arbitrary  random  phase.  Equation  (C.4)  may  also  be 
o 

written  as  follows: 
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1 


v  -  v  +  v  (u  -  u  )  +  (v  -  V  )(u  -  u  )  +  ..  (v  -  v  )(u  -  u  ) 

n  o  o  n+1  o  1  o  n+1  1  n  n-1  n+1  n 


v  +v(u  -u)+^(v-v  )(u  -  u  ) 

o  o  n+1  o  ^ —  1  i-1  n+1  1 

i-1 


(C.5) 


Each  term  in  the  above  expression  is  independent  and  of  zero  mean 
with  boundary  conditions  specified  in  equations  (C.2b)  and  (C.2c). 
Hence  the  variance  may  be  written: 


var{v  } 


2  2 
c  or  +  c  o"  (u 
32  2  2  n+1 


2  sr  2 

-  u  )  +  >  or  (u  -  u 
o  2  i  i-1 


)(u  -  u  ) 

n+1  i 


(C.6) 


We  can  evaluate  the  above  summation  in  the  limit. 


n  ->  oo  £— 

max{u  -u  }  ->  0  i-1 

i  i-1 


n  n 

Z  2  2  2  ( 

<T  (u  -  u  )(u  -  u  )  -  <r  V 

2  i  i-1  n+1  i  2  ) 


n+1 

2  2  (  2 
)  -  <r  \  (u  -  y)  dy 
2  )  n+1 


(u  -  u  ) 
2  n+1  o 


Denoting  points  in  time  as 


u  -  u  and  u  -  u  +t 

o  n+1 


the  limit  of  equation  (C.6)  can  be  written 
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3 


Var {v(u+t)} 


-  I  I  |V|  Var{v  } 
n  — >  <xd  n 

max{u  -u  }  ->  0 
i  i-1 


2  2  t 

O'  (c  +  c  t  +  —  ) 
2  3  2  3 


2  2  2  2  t 

ccr  +  co't+cr  — 
3  2  2  2  2  3 


(C.7a) 


Ue  can  also  evaluate  the  variance  for  t«0 


Var{v(u)}  *  c  cr 
3  2 


(C.7b) 


From  the  Independence  of  the  values,  v  and  v  ,  and  the  Independence 

o  o 

of  the  Increments,  we  can  evaluate  the  covariance 


u+t 

*  (  .  *  * 

{v(u+t)v  (u)}  -  E{ (v  +  \  v(y)dy)(v  )}  ■  E{v  v  }  -  var{v(u)} 

0  f  o  o  o 


cov 


2  * 

c  <T  -  cov{v  (u+t)  v(u)} 
3  2 


<C.7c) 


Equation  (C.7)  provides  the  elements  of  the  covariance  matrix,  H  , 

v 

for  the  non~statlonary  process 


v(u+t) 

v(u) 


The  covariance  matrix  Is  obtained  by  making  the  appropriate 
substitutions  Into: 
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M 

v 


v(u+t) 

v(u) 


[ 


* 

V  (u+t) 


* 

var{v(u+t)}  cov{v(u+t)v  (u)} 

* 

cov{v  (u+t)v(u)}  var{v(u)} 


(C.8) 


Substitution  of  equation  (C.7)  into  (C.8)  resulta  in  the  desired 
expression  which  is  used  in  equation  (2.14)  of  the  text.. 

C.2 .  Integration  Formula  Approach 

Next  we  shall  use  an  integration  formula  for  determining  the 
covariance  matrix  of  the  phase  random  process  resulting  from  the  random 
walk  frequency  process.  We  first  determine  the  variance  of  the 
integration  process,  equation  (C.l),  at  some  Instant  of  time,  say  (u+t). 

From  [C.4]  (equations  (5.3.1)  and  (5.3.4)),  equation  (C.l),  the 
condition  given  in  equation  (C.2b)  and  (C.2c),  and  the  independence 
assumption  given  above,  we  have  for  the  integral  of  the  random  walk 
frequency  process 


I. 

v 


u+t 


5 


v(y)dy 


and  for  the  variance 
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2 

E{| I. I  } 
v 


2 

<r 

2 


u+t 


c  t  dq 
2 


+ 


u+t 


u 


-  q 

^  (p-u)  dp 

—  u 


u+t 


q 


2 

r 

2 


to 

q  2 

(p-u) 

u+t  - 

c  t  +  \  dq 

2  / 

-  + 

2 

(q-u)p 

u 

u 

q  -1 

u+t 

r~ 

2 

“71 

2 

c  t  + 
2 


dq 


u 


(q-u) 


+  (q-u)(u+t-q) 


3 

2  2  t 

<T  (c  t  + - ) 

2  2  3 


Subs tit ion  of  equation  (C.ll)  into  equation  (C.9)  results  in 


3 

2  2  t 

var{v(u+t)}  «<r(c+ct  +  —  ) 
2  3  2  3 


u)  dp 


> 


(C.ll) 


(C. 12a) 
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This  is  the  s*-  expression,  equation  (C.7a),  obtained  using  the 
approximation  function  in  approach  one  above.  It  can  be  shown,  using 
the  same  reasoning  as  presented  in  approach  one  that  equation  (C.7b)  and 
(C.7c)  still  hold. 


2 

var{v(u)}  -  c  o 
3  2 

*  2 
cov{v(u+t)v  (u)}  *  C  O' 

3  2 


(C.12b) 


(C.12c) 


C.3 •  State  Space  Approach 

We  now  use  a  state  variable  approach  for  determining  the  variance, 
var{v(u+t)} . 

The  procedure  we  follow  is  outlined  in  [C.5]  (section  6.3.1, 
Property  14,  pp  533).  It  states: 

The  variance  matrix  of  the  state  vector,  v(t)  of  a  system 

v(t)  =■  F(t)  v(t)  +  G(t)  u(t)  (C.13a) 

satisfies  the  differential  equation 

T  T 

A  (t)  -  F(t)A  (t)  +  A  <t)F  (t)  +  G(t)Q(t)G  (t)  ( C •  1 3b ) 

V  V  V 

with  condition 
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A  (u)  -  E[  v(u)  v  (u)  J 


(C.13c) 


Observe  that  A  ( t )  «  K  (t,t)  is  the  covariance  matrix  of 


the  vector  v. 


To  use  equation  (C.13)  for  the  random  walk  frequency  we  set 


v  (t) 
1 

v  (t) 
2 


0  1 


0  0 


A  (t)  -  E[  v ( t )  V  (t >] 


where  v  (u)  denotes  the  phase  random  process,  v  (u)  denotes  the  frequency 
1  2 
random  process  and  u(t)  is  a  zero  mean  white  noise  process  with 

variance,  Q(t)  ■  1.  This  latter  result,  Q(t)  is  obtained  from  [C.5] 

(Property  13,  equation  (257)  and  the  result  immediately  following 

equation  (254a)). 

Substitution  of  the  above  into  equation  (C.13a)  and  (C.13b)  gives 
us  (C.14a)  and  (C.14b)  respectively 


v  (t)  v  (t) 

1  2 


v  (t) 
2 


<T  u(t) 
2 


(C.l4a) 
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r—  — 

A  - 

V 

-l  -1— 

V  V  V 

1  1  2 

•  • 

V  V  v1 

12  2 

- 

0  1 

0  0 

r““ 

v  a  V  V 

1  12 

V  V  v  1 

12  2 

+ 

V  *  V  V 

1  12 

-  —z 

V  V  V 

12  2 

0  0 

1  0 

r-"  _ 

i—  — 

2  v  v  v  z 

0  0 

12  2 

+ 

v  1  0 

2 

0  o' 

2  J 

—  7  - 

(C.14b) 


where  the  complex  notation,  *,  has  been  omitted  because  the  process, 

v(t).  Is  a  phase  random  process  which  has  been  modeled  as  a  real  process 

throughout  this  dissertation.  Dependence  on  time,  v  -  v(t).  Is  Implied, 

and  notation,  E[v]  ■  v  ,  is  adopted,  v  denotes  the  phase  random 

1 

process  and  v  denotes  the  frequency  random  process,  with  conditions 
2 

at  time,  u,  taken  from  equations  (C.2c)  and  (C.2b)  respectively  In 
conjunction  with  the  independence  relation  between  v(u)  and  v(u).  The 
conditions  at  time,  u,  equation  (C.13c),  may  be  stated 


T 

A  (u)  -  E[  v(u)v  (u) 
v 


(C.l4c) 
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The  variance  matrix  A  is  easily  solved.  From  equation  (C.14)  we  can 
collect  three  first  order  differential  equations  along  with  their 
boundary  conditions. 


Differential  Equation  Boundary  Condition 


v^«2vv 
1  12 


vi(u)  ■  c  <r 
1  3  2 


an  arbitrary  constant 


vv  ■»  v  v  ■  v1 
12  12  2 


V  V  (u) 
1  2 


-  0 


vz  -  a- 
2  2 


v\u) 

2 


2  a  function  of  u, 

c  <r  :  long  term  drift 
2  2  offset  frequency 


Solving  for  v^(t)  first,  we  have 
2 


’ 


vZ(t)  dt  +  k 


$ 


2  2 
or  dt  +  k  -  cr  t+k 
2  12  1 


Applying  its  boundary  condition  at  time,  u,  we  can  solve  for  k  , 

1 


—l 


2  2  2 
v*(u)  -  <r  u  +  k  -  c  O'  — k  -  (c  -  u)cr 
2  2  122  122 


Substituting  for  k  ,  we  obtain  the  solution  for  v  (t), 
1  2 


—  2 

v  (t)  -  <r  (t  +  c  -  u) 
2  2  2 


Next,  we  solve  for  the  quantity,  v  v  (t). 

1  2 
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r  v  (t)  -  \  v 

12  )  1 


v  (t)  dt  +  k 
2  2 


5 


_  2  t 

v*(t)  dt  +  k  -  O'  (  —  +  [c  -  u]  t)  +  k 
2  2  2  2  2  2 


Applying  Its  boundary  condition  at  time,  u,  we  can  solve  for  k  , 

2 


2  2 

_  2  u  2  u 

vv(u)«<r  (cu - )  +  k  -  0  — »  k  -  cr  ( - cu) 


1  2 


2  2  2 


2  2  2 


Substituting  for  k  ,  we  obtain  the  solution  for  v  v  (t) 
2  12 


_  2  [t-u] 

v  v  (t)  -  <r  (  - +  c  [t-u]  ) 

1  2  2  2  2 


Finally,  we  solve  for  the  quantity,  v2(t) 

1 


r^(t)  -  f  v2i 

1  )  1 

r*(t)  -  2  )  ~ 
1  /  1 


<t)  dt  +  k 


v  dt  +  k 
2  3 


3  2 

2  [t-u]  c  [t-u] 

2cr  ( - + - )  +  k 

2  6  2  3 


Applying  the  boundary  condition  at  time,  u,  we  can  solve  for  k 


-2,  2  2 

v  vu )  -  k  -  c  cr  — ^  k  -  c  cr 

1  3  3  2  3  3  2 


Substituting  for  k  ,  we  obtain  the  solution  for  v*(t) 
3  1 
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1 


I- 

i 

I. 


3 

_  2  [t-u]  2 

v  (t)  ■  a  ( - +  c  [t-u]  +  c  ) 

1  2  3  2  3 


The  variance  of  the  phase  random  process,  v  (t),  is  obtained  from  this 

1 

last  expression.  After  a  change  in  variables,  t'  ■  (t-u),  we  have 

3 

2  _  2  2  t' 

var{v  (u+t')}  ■  E{v  (u+t')}  ■  v2(u+t')  -  <r  (c  +  c  t'  + - )  (C.15a) 

1  1  1  2  3  2  3 


This  is  the  same  expression  obtained  in  equation  (C.7a)  of  approach 
one  and  equation  (C.12a)  of  approach  two. 

Using  the  lndepencence  assumption  we  have 


2 

var{v(u)}  -  c  <T 
3  2 

*  2 

cov{v(u+t)v  (u)}  ■  c  <T 

3  2 

which  are  the  desired  results. 


(C.15b) 


(C.15c) 
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C . 4 .  Summary  of  Approaches 


From 

function, 

frequency 


approach  one,  two,  or  three  above,  we  obtain  the  covariance 

M  of  the  phase  random  process,  v(u),  resulting  from  the 
v 

random  walk  process,  v(u); 


*T 

M  =  E{  v  v  } 


{ 


v(u+t) 

v(u) 


c 


*  * 

V  (u+t)  V  (u) 


2 

=  or 
2 


3 

2  t 

c  +  c  t  +  — 
3  2  3 


c 

3 


c 

3 


(C. 16) 


This  result  is  used  in  the  text,  cf.  section  II. 
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APPENDIX  D 

Determination  of  Second  and  Fourth  Order  Moments 
of  the  RF  Signal  Process 

0.0.  Introduction 

There  are  perhaps  many  approaches  that  can  be  used  for  determining 
the  higher  order  moments  of  oscillator  rf  signal  process  due  to  white 
phase,  y(u),  random  walk  phase,  x(u),  and  random  walk  frequency,  v(u). 
In  this  appendix,  we  use  two  approaches,  the  characteristic  function  and 
a  method  involving  independent  increments. 

D.l .  Characteristic  Function 
D.1.1.  Near  White  Phase,  y(u) 

We  determine  the  higher  order  moments  of  b(u):  t  »  s,  for  y(u) 
Gaussian  distributed  with  the  aid  of  figure  D.l  and  the  characteristic 
function  method  [D.l]  (chapter  4  and  8),  [ D - 2 ]  presented  in  Appendix  A. 

We  define/construct  the  four  element  vector 

T 

£  [  y(u+s+t)  y(u+t)  y(u+s)  y(u)  ] 

4 

Using  figure  D.l  we  can  easily  write  the  covariance  matrix  for 
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• 

c 

c  -  2s/T  c  -  1 

c  -  1 

T 

o 

o  o  o 

o 

o 

:  s  <  — 

c  -  2s/T 

c  c  -  1 

c  -  1 

2 

r 

2 

H  -  <r 

0  o  o  o 

o 

^4  o 

c  -  1 

c  -  1  c 

c  -  2s/T 

T 

o 

o  c 

o  o 

o 

:  t  >  —  +  s 

c  -  1 

c  -  1  c  -  2s/T  c 

2 

o 

o  o 

o  o 

J 

(D.l) 

r 

— 

c 

c  -  1  c  -  1 

c  -  1 

T 

o 

o  0 

0 

o 

- 

c  -  1 

c  c  -  1 

c  -  1 

:  s  >  — 

2 

2 

o 

o  o 

o 

■ 

M  -  <r 

£4  o 

c  -  1 

c  -  1  c 

c  -  1 

T 

o 

O  0 

o 

o 

:  t  >  —  +  s 

'W 

c  -  1 

c  -  1  c  -  1 

c 

2 

o 

o  o 

0 

We  will  denote  the 

elements  of  the 

matrix,  M 

by  m  .  Substitution  of 

y4  lj 

equation  (D.l)  into  equation  (A. 3) 

leads  to  the  fourth  order  moment 

■HBi 

for  the  rf  waveform,  b  (u). 

I 

y 

2 

R 

(t)  - 

exp(-  4s<r  /T  ) 

:  s 

<  T  /2  <  t-s 

*4>y 

o  0 

o 

2 

(D.2a) 

'  (■ 

R 

(t)  - 

exp(-  2cr  ) 

:  T 

/2  <  s  <  t-T  /2 

j>n 

0 

0 

o 

In  much  the  same  manner  the  second 

order  moment  of  the  rf  waveform  is 

• 

2 

R  (s)  - 

exp(-2scr  /T  ) 

:  s 

<  T  /2  <  t-s 

y 

o  o 

0 

(D.2b) 

- 

2 

is 

R  (s)  - 

exp  (-o’  ) 

:  T 

/2  <  s  <  t-T  /2 

y 

o 

o  o 

. 

•% 

150  - 

-  150  - 


1 


The  fourth  and  second  order  moments  for  the  white  phase  component 


are  given  by  equations  (D.2a)  and  (D.2b)  respectively. 


D.1.2.  Random  Walk  Phase,  x(u) 


We  revisit  the  random  walk  phase  process.  Our  vector  random 


process  is 


x  ^  [  x(u+s+t)  x(u+t)  x(u+s)  x(u)  ] 
4 


With  the  aid  of  figure  D.2  we  write  the  covariance  matrix 


c+s+t  c  +  t  c+s 

1  1  1 

c+t  c  +  t  c  +  t 

2  11  1 

r 

1  c  +  s  c  +  t  c+s 
1  1  1 


:  lt|  <  I s | 


c  +  s  +  t  c  +  t  c  +  s 

1  1  1 

c  +  t  c  +  t  c  +  s 

2  1  1  1 
M  =  rr 

x4  1  c  +  s  c  +  s  c  +  s 

111 


(D.3) 


Is  I  <  lt| 


Substitution  of  equation  (D.3)  into  equation  (A. 3)  leads  to  the 
fourth  order  moment  of  the  rf  signal  when  the  phase  is  a  random  walk 
Gaussian  process. 
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FIGURE  D.l  AUTO-CORRELATION  FUNCTION  OF  NEAR  WHITE  PHASE,  y(u) 

SLOPE*  0*,2 


FIGURE  D.2  VARIANCE  OF  RANDOM  WALK  PROCESS,  x(u) 
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2 

R  (t)  -  exp(-<f  1 1 1  ) 

tt*  1 

2 

R  (t)  “  exp  (-o’  J  s  |  ) 

H*  1 


Itl  <  is! 


Isl  <  Itl 


(D.4a) 


This  is  in  agreement  with  equation  (D.9)  which  will  be  derived  as 
the  product  of  two  independent  increments.  The  second  order  moment  is 
easily  obtained  using  the  same  procedure, 


2  |s| 

R  (s)  ■  exp(-<r - )  (D.4b) 

x  1  2 


D.1.3.  Random  Walk  Frequency,  v(u) 

Four  Element  Vector 

In  order  to  determine  the  fourth  order  moment  of  b(u)  when  the 
phase  random  process  arises  from  a  Gaussian  random  walk  frequency  com- 
ponent,  v(u),  we  must  determine  the  covariance  matrix  of 

T 

v  -  [  v(u+s+t)  v(u+t)  v(u+8)  v(u)  ] 

4 

as  we  did  for  the  white  phase  and  random  walk  phase  processes  above. 

The  approximating  function,  [D.3],  (chapter  7,  sec. 5),  a  discrete 

approximation  of  the  continuous  function,  can  be  made  to  provide  a 

systematic  procedure  for  determining  the  elements  of  the  covariance 

matrix,  M  .  Using  figure  D.2  for  the  case,  |s|  <  |t| 
v4 
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>4 


m 


v 

m 


v+v(u-u)  + 
o  o  m  o 


i-1 


(  v  -  v  )(  u  -  u  ) 
i  i-1  m  i-1 


(D.5) 


where  m  -  o,  n,  N,  n+N  is  the  number  of  discrete  steps  constituting 

time  u  -  u  ,  u+s,  u+t,  u+t+s  respectively,  see  figure  D.2. 
o  o  o  o 


For  v  ,  v  ,  and  (  v  -  v  )  zero  mean  independent  and  real,  we  may  write 
o  o  i  i-1 

the  covariance  matrix  elements  of  v  . 

“4 

The  diagonal  terms  can  be  evaluated  as  follows. 

2 

Var{v  }  A  E{  v  } 
m  m 

Substitution  of  equation  (D.5)  intq  the  above  we  have 


m 

2  .2  2  ^ —  .  .  2  2 

Var{v  }  -  E{v  }  +  E{v  (u-u)}+\  E{(v-v  )(u-u  )} 

m  o  o  m  o  ^ —  i  i-1  m  i-1 

i-1 


Evaluating  the  expectation  and  using  the  property  of  the  drift  variance 
and  stationary  independent  increments  associated  with  the  Brownian 
motion  process  under  the  summation  we  have 


Var{v  } 
m 


m 

2  2  2  2  2 
c<r  +  cff(u-u)+^  <T(u-u  )(u-u  ) 

3  2  2  2  m  o  2  m  i-1  i  i-1 

i-1 


which  in  the  limit  leads  to 
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Var{v(u  +t)} 
o 


I  M 

o  <  m  — ^  qcd 

max{u  -  u  }  ->  o 
i  i-1 


Var{v  } 

m 


u  +t 

2  2  2  o/2  2 

c  <r+  c  <r  t  +  1  <r  (u  +  t  -  x)  dx 

3  2  2  2  )  2  o 


2  2  t 

cr  (c  +  c  t  +  —  ) 
2  3  2  3 


(D.6a) 


where  x  is  the  variable  of  intergration  corresponding  to  the  dummy 

index,  i,  and  u  +  t  is  designated  as  time,  u  . 

o  it 

The  off-diagonal  terms  are  similarly  evaluated.  Let 

m  <  m  denote  time,  u  <  u  :  where  ml,  m2  are  the  number  of  discrete 
12  ml  m2 

Increments  along  the  time  axis,  u,  and  u  ,  u  correspond  to  shifts 

ml  m2 

from  u  along  u. 
o 


Next  we  evaluate  the  off  diagonal  elements  of  the  covariance  matrix. 


Cov{v  v  }  ^  E{v  v  } 
ml  m2  ml  m2 

Substitution  of  equation  (D.5)  into  the  above  gives 


2  .2 

Cov{v  v  }  -  E{v  }  +  E{v  (u  +  u  -  u  )(u  +  u  -  u  )} 
ml  m2  o  o  o  ml  o  o  m2  o 


where  ml,  and  m2  denotes  any  of  the  subscripts,  m  -  o,  n,  N,  n+N. 


-  155  - 


From  the  stationary  independent  increments 


Cov  {v  v  } 
ml  m2 


ml 


2  2 

cff+cffu  u  + 
3  2  2  2  ml  m2 


X 


2 

E{(v  -  v  )  (u  - 
i  i-1  ml 


u  )(u  - 

i-1  m2 


u  )} 

i-1 


ml 

2  V  2 

-  <r  (c  +  C  U  U  )  +  >  c  (u  -u  )(u  -  u  )(u  -  u  ) 

2  3  2  ml  m2  2  ml  i-1  m2  i-1  i  i-1 

i-1 


Taking  the  limit  as  the  number  of  discrete  steps, 
approach  infinity  we  have 


for  m  and  m  ,  each 
1  2 


Cov{v(u  +  u  )v(u  +  u  )} 
o  ml  o  m2 


ml  <  m2  ->  od 

max {u  -  u  }  ->  0 
i  i-1 


Cov {v  v  } 
ml  m2 


2 


u  +  u 
o  ^  ml 


«  O'  (c  +  c  u  u  )  + 
2  3  2  ml  m2 


<r  (u  +  u  -  x)(u  +  u  -  x)  dx 
2  o  ml  o  m2 


u 


o 


2  3 

u  u 

2  ml  ml 

-c(c+cu  u  +  u  - —  -  - ) 

2  3  2  ml  m2  m2  2  6 


(D.6b) 


Equation  (D.6)  can  be  used  to  determine  the  covariance  matrix  of  the 

vector,  v  ,  equation  (D.6a)  for  the  diagonal  elements  and 
4 

equation  (D.6b)  for  the  off-diagonal  elements.  For  example, 

substituting  t  for  u  and  t+s  for  u  we  can  determine  the  covariance 

ml  m2 
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matrix  element,  m  ;  or  substituting  for  u  and  t  for  u  we 
12  ml  m2 

determine  the  element  m  .  Using  the  vector  definition  from  above 

23 


T 

V  ■  [  v(u+s+t)  v(u+t)  v(u+s)  v(u)  ] 
”4 


the  covariance  matrix  element,  m  can  be  determined  from  equation  (D.6b) 

23 


m 

23 


Cov{v(u  +  u  )v(u  +  u  )} 
o  ml  o  m2 


U  «  8,  U  -  t 

ml  m2 


3 

u 

ml 


|U  «  8,  U  -t 

ml  m2 


2  3 

2  s  s 

■  <r  (c  +  c  st  + —  t - ) 

2  3  2  2  b 


2 

u 

2  ml 

■<y(c+cu  u  +  u  — — 

2  3  2  ml  m2  m2  2 


Substitution  of  the  appropriate  variables  for  u  and  u  into  equation 

ml  m2 

(D.6)  for  the  four  element  vector,  the  covariance  matrix  elements  are 
determined  as  a  function  of  s  and  t.  Then  substituting  these  matrix 
elements  into  the  characteristic  function  of  equation  (A. 3a),  the  fourth 
order  moment  of  the  rf  waveform,  b(u),  is  found. 


2 

2  s 

R  (t)  «  expC-o" - f3t  -  s]) 


:  Isj  <  |t|  (D.7a) 


rr 


N 


In  a  similar  fashion  the  second  order  moment  of  b(u)  is  determined  using 
equation  (A. 3b). 

2 

2  s 

R  (s)  ■  exp(-o"  —  [3c  +  s])  (D.7b) 

b  2  6  2 

Next  we  use  an  alternative  approach  for  determining  the  phase  covariance 
matrix. 

D.2.  Independent  Increments 

This  approach  makes  use  of  the  independent  increments  which 
constitute  the  random  walk  processes  Involving  phase  and  frequency.  It 
has  the  advantage  of  decomposing  a  four  by  four  element  matrix  and 
corresponding  four  element  vector  characteristic  function  into  a  two  by 
two  matrix  and  two  element  vector  characteristic  function.  We  will  find 
that  the  procedure  greatly  simplifies  determination  of  the  covariance 
matrix  for  the  random  walk  phase  and  frequency  processes  as  well  as 
serving  as  a  check  against  the  method  for  the  four  element 
characteristic  function  just  presented. 

D.2.1.  Stationary  Increments  (Brownian  Motion/ Random  Walk  Phase) 

For  any  random  walk  phase  process  the  fourth  order  moment  required 
in  the  theorems  of  chapter  III  can  be  evaluated  as  follows. 

*  * 

R  (t)  -  E[  b(u+s+t)  b  (u+t)  b  (u+s)  b(u)  ]  (D«8) 

>*  E[  exp{  jx(u+s+t)  -  jx(u+t)  -  jx(u+s)  +  jx(u)  }  ] 
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By  examining  the  exponent  and  making  use  of  the  stationary  independent 
increments  the  phase  term,  x(u),  in  the  exponent  can  be  rewritten 
as  the  sum  of  two  independent  random  variables. 


[x(u+s+t)  -  x(u+t ) ]  +  { x (u )  -  x(u+8)]  :  |s|  <  |t| 

x(u+s+t)  -  x(u+s)]  +  [x(u)  -  x(u+t)]  :  I s |  >  |t| 

This  is  illustrated  in  figure  D.3. 


The  fourth  order  moment  for  the  random  walk  Gaussian  distributed 
phase  component  only  is 


2  2 

tr  cr 

112 

R  (t)  -  exp( - |t|)  exp( - |t|)  «  exp(-cr  |t|)  |t|  <  |s| 

0b  2  2  1 

(D.9) 

2  2 

O  O' 

1  1  2 

R  (t)  **  exp( - I s |  )  exp( - |s|)  “  exp(-cr  |s|)  :  I s I  <  lt| 

0b  2.2  1 


where  |s|,  |t|  is  the  appropriate  first  order  increment  which  emerged 
from  Appendix  B.  The  result  is  identical  to  the  fourth  order  moment 
given  in  equation  (D.4a),  where  the  characteristic  function  method  was 
used. 


D.2.2.  Two  Element  Increment  Vector  for  Random  Walk  Frequency 


Let  us  refer  back  to  figure  D.2.  We  define  the  two  element  vector 


T 

v  ^  [v(u+s+t>  -  v(u+t)  v(u+s)  -  v(u)] 
2  “ 
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FOR  CORRELATION  FUNCTIONS  HAVING  FACTORS  E.O.  (EXPONENTIAL  ORDERK-O 
e.g.  RANDOM  WALK  PHASE  &  RANDOM  WALK  FREQUENCY 
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0  0 
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U0  4>+s  ty-t  +  s 

R^  (t):«"o,12  |S( 


FIGURE  D.3  DOMAIN  FOR  FOURTH  ORDER  MOMENTS  OF  WAVEFORM,  b(u) 


Again,  in  our  approximating  function  approach,  we  set  a  number  of 

discrete  increments  to  ml  and  m2  where  ml  <  m2  and  substitute 

u  for  the  shift,  t  and  u  for  the  shift,  t+s,  and  u  for  the 
ml  m2  m2-ml 

shift,  s.  The  shift  is  from  u  along  the  time  axis,  u.  The  diagonal 

o 

terms  of  the  covariance  matrix  of  the  vector  v  can  be  written 

~2 


Var{v  -  v  }  ^  E{(v  -  v  )  } 
m2  ml  m2  ml 


Substitution  of  the  approximating  model,  equation  (D.5),  into  the  above 
gives  us,  after  some  manipulation, 


•  2  2 
Var{v  -  v  }  *  E{v  (u  -  u  )  } 
m2  ml  o  m2  ml 


ml 


m2 


r—  'c* —  .  .  — »2  jr—  v —  .  .  — 

+  E{  >  (V  -  V  )(U  -  U  )  +  >  (V  -  V  )(u  -  U  ) 

L_  /—  i  i-1  m2  ml— J  L-  i  i-1  m2  i-1  -J 


i-1 


i*ml+l 


Taking  the  limit  we  have 


Var{v(u  +  u  )  -  v(u  +  u  )}  - 

m2  ml 


I  M 

ml,  m2  — >  a® 


Var{v  -  v  } 
m2  ml 


max{u  -  u  }  ->  0  *— 
i  i-1 


u  +  u 
o  ml 


2  2 
c  <r  (u  -  u  )  + 
2  2  m2  ml 


u  +  u 
o  m2 


2  2 
<r  (u  -  u  )  dx  + 
2  m2  ml 


2  2 
cr  (u  +  u  -  x)  dx 
2  o  m2 


u 


u  +  u 
o  ml 


The  desired  result  is 


Var{  v(u  +  u  )  -  v(u  +  u  )} 
m2  ml 


(u  -  u  )  -» 

2  m2  ml  | 


CQU  -U  )  +  U  (U  “U 

2  m2  ml  ml  m2  ml 


(D.lOa) 


The  off  diagonal  terms  may  be  found  by  designating  the  second 

element  in  the  vector  v  by  the  substitution,  u  for  u  +  a  where 

2  m2-ml  o 


0<u  -u  <u  <u 

m2  ml  ml  m2 


Then  it  follows 


Cov{(v  -  v  )(v  -  v  )}  ^  E{(v  -  v  )(v  -  v  )} 

m2  ml  m2-ml  o  m2  ml  m2-ml  o 


m2-ml 

.  2  r —  .  •  2 

E{v  (u  -  u  )(u  -  u  )}  +  E  \  (v-v  )  (u  -u  )(u  -  u  ) 

o  m2  ml  m2-ml  o  ^ —  i  i-1  m2  ml  m2-ml  o 

i-1 


Taking  the  limit  we  have 
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Cov{[v(u  +u  )  -  v(u  +  u  )][v(u  +  u  )  -  v(u)]} 
m2  ml  m2-ml 

I  i  M  2 
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u  )(u  +  u  -  x)  dx 
ml  o  m2-ml 


(D.lOb) 


Making  the  appropriate  substitutions  for  u  and  u  Into  equation 

ml  m2 

(D.10)  for  the  two  element  Increment  vector,  the  covariance  matrix 
elements  are  determined  as  a  function  of  time  shifts  s  and  t. 
Substitution  of  equation  (D.10)  Into  equation  (A. 3b)  results  in  the 
fourth  order  moment  of  the  rf  signal,  b(u).  Note  that  it  Is  identical 
to  a  previous  result  found  in  equation  (D.7a), 


2 

R  (t)  -  expP-  6 - (3t  -  s)l  s  I  s  |  <  1 1 1  (D.7a ) 

oob  w-  2  6 


The  second  order  moment  Is  determined  by  similar  methods  using  a  one 
dimensional  vector  for  v(u),  i.e.  a  scaler 


R  (s) 


2 

- (3c  +  s)l 


(D.7b) 


D.3.  Summary 


Using  characteristic  functions  a  systematic  procedure  was  obtained 
[D.2]  for  evaluating  moments  and  auto-correlation  functions  of  the  rf 
signal,  b(u),  as  required  in  the  theorems  of  chapter  III.  This  was 
demonstrated  for  each  phase  instability  process  treated  here;  the  white 
and  random  walk  phase  and  frequency  including  initial  starting  phase. 
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